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The non-equilibrium dynamics of classical random Ising spin chains with non-conserved magne- 
tization are studied using an asymptotically exact real space renormalization group (RSRG). We 
focus on random field Ising spin chains (RFIM) with and without a uniform applied field, as well as 
on Ising spin glass chains (SG) in an applied field. For the RFIM we consider the universal regime 
where the random field and the temperature are both much smaller than the exchange coupling. In 
that regime, the Imry-Ma length that sets the scale of the equilibrium correlations is large and the 
coarsening of domains from random initial conditions (e.g. a quench from high temperature) occurs 
over a wide range of length scales. The two types of domain walls that occur diffuse in opposite 
random potentials, of the form studied by Sinai, and domain walls annihilate when they meet. Using 
the RSRG we compute many universal asymptotic properties of both the non-equilibrium dynamics 
and the equilibrium limit. We find that the configurations of the domain walls converge rapidly 
towards a set of system-specific time-dependent postitions that are independent of the initial condi- 
tions. Thus the behavior of this non-equilibrium system is pseudodeterministic at long times because 
of the broad distributions of barriers that occur on the long length scales involved. Specifically, we 
obtain the time dependence of the energy, the magnetization and the distribution of domain sizes 
(found to be statistically independent). The equilibrium limits agree with known exact results. We 
obtain the exact scaling form of the two-point equal time correlation function (So(t)S x [t)} and the 
two-time autocorrelations (So(t')So(t)). We also compute the persistence properties of a single spin, 
of local magnetization, and of domains. The analogous quantities for the ± J Ising spin glass in an 
applied field are obtained from the RFIM via a gauge transformation. In addition to these we com- 
pute the two-point two-time correlation function {So(t)S x (t)){So(t')S x (t')) which can in principle be 
measured by experiments on spin-glass like systems. The thermal fluctuations are studied and found 
to be dominated by rare events; in particular all moments of truncated equal time correlations are 
computed. Physical properties which are typically measured in aging experiments are also studied, 
focussing on the response to a small magnetic field which is applied after waiting for the system to 
equilibrate for a time t w . The non-equilibrium fluctuation-dissipation ratio X(t,t w ) is computed. 
We find that for (t — t w ) ~ t", with a < 1, it is equal to its equilibrium value X = 1 although time 
translational invariance does not hold in this regime. It exhibits for t — t w ~ t w an aging regime 
with a non-trivial X = X(t/t w ) 1, but the behaviour is markedly different from mean field theory. 
Finally the distribution of the total magnetization and of the number of domains is computed for 
large finite size systems. General issues about convergence towards equilibrium and the possibilities 
of weakly history-dependent evolution in other random systems are discussed. 
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I. INTRODUCTION 



In many systems, the development of long range order is controlled by the dynamics of domain walls. The coarsening 
of domain structures evolving towards equilibrium has been studied extensively in pure systems [1,2] but little is known 
quantitatively about domain growth in the presence of quenched disorder. In random systems, the nonequilibrium 
dynamics plays an even more important role than in pure systems and is most relevant for understanding experiments, 
since many random systems become glassy at low temperatures with ultraslow dynamics which prevent full thermal 
equilibrium from being established within the accessible time scales. This slow dynamics is due, at least in cases 
in which it is qualitatively understood, to the large free energy barriers that must be overcome for order to be 
established on long length scales. Some of the best studied cases both experimentally and theoretically are a variety 
of random magnetic systems, particularly spin glasses and random field magnets [3-5]. Both of these systems have 
engendered a great deal of controversy about their equilibrium behavior and the resolution of these controversies 
has been greatly hampered by the inability of experiments to reach equilibrium. One might well argue, however, 
that the most interesting properties of such random systems are, in fact, not their equilibrium properties, but the 
non-equilibrium dynamics involved in their Sisyphean struggle to reach equilibrium. 

Because of the dominance of the behavior of so many random systems by the interplay between equilibrium and 
non-equilibrium dynamic effects, it is important to find models with quenched randomness for which solid results about 
dynamics can be obtained. In particular, one would like to be able to analyze the effects of activated dynamics caused 
by large barriers and compare results to predictions which come from either phenomenolgical scaling approaches to 
the dynamics — often known somewhat misleadingly as "droplet models" [5] , or from mean field approaches [6] . Many 
of the interesting phenomena go under the general name of aging — the dependence of measurable properties such as 
correlations and responses on the history of the system, in particular on how long it has equilibrated for: its "age" . 
For example, one would like to compute quantities which probe the violations of the fluctuation dissipation relations 
caused by non-equilibrium effects and compare them with results obtained in mean field models [6]; this has been 
done previously for coarsening of pure models [7-9] . 

One of the simplest non-trivial random models, and one in which coarsening occurs, is the random field Ising model 
(RFIM) in one dimension (ID). Although this system does not have true long-range order, for weak randomness it 
exhibits a wide range of scales over which dynamics is qualitatively like those of other random systems, especially those, 
like two dimensional random field magnets and two-dimensional spin glasses, which do not have phase transitions 
but exhibit much dynamic behavior qualitatively similar to their three-dimensional counterparts which do have phase 
transitions. In particular, the weakly random ID RFIM has a wide range of length scales over which the typical size of 
ordered domains grows logarithmically with time at low temperatures. In ID, the RFIM is equivalent to a spin glass 
in an applied magnetic field; this, or some analogous ID systems should be conducive to experimental investigation. 
For a recent review on the RFIM see e.g. [11]. 

The equilibrium properties of the ID RFIM and of the ID spin glass in a field have been extensively studied [11-25]. 
Several thermodynamic quantities such as the energy, entropy and magnetization have been computed exactly at low 
temperature for a binary distribution of the randomness in [14-16] and for continuous distributions in [23-25]. Results 
are also available for distributions of bonds with anomalous weight near the origin [17]. The free energy distribution 
has been studied in [18]. Equilibrium correlation functions are harder to obtain and only a few explicit exact results 
exist. For the binary distribution some results are presented in [21]. Certain special limits have been solved, such 
as the infinite field strength limit which is related to percolation [25], but for, e.g., higher cumulants of averages of 
truncated correlations, only the general structure has been discussed [22]. 

In this paper we derive a host of exact results for the nonequilibrium dynamics of the RFIM and for the ID Ising 
spin glass in a field, and, as a side benefit, also obtain many equilibrium results, some new. Those of our results for 
equilibrium quantities that can be compared with previously known ones are found to agree with these. 

Here we focus on the universal regime of the RFIM in which both the random field and the temperature are much 
smaller than the exchange coupling, the length and time scales are sufficiently long and the random field dominates 
the dynamics. With the random field being weak and the temperature low, the equilbrium correlation length is long 
and essentially the same as that at zero temperature — the Imry-Ma length. The coarsening of domains starting 
from initial conditions with only short distance correlations — such as a quench from high temperature — will, after 
a rapid initial transient, be dominated by the randomness over a wide range of time scales . 

The basic tool that we use in this paper is a real space renormalization group (RSRG) method which we have 
developed recently to obtain exact results for the non equilibrium dynamics of several ID disordered systems [26]. 
Most of our previous results have been for the Sinai model that describes the diffusion of a random walker in a ID 
random static force field which is equivalent to a random potential that itself has the statistics of a ID random walk 
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[27]. It can readily be seen that individual domain walls in the classical RFIM diffuse in a random potential of exactly 
this Sinai form, the complication being that they annihilate upon meeting. As shown in [26,29] the RSRG can also 
be applied to many-domain-wall problems such as that which corresponds to the RFIM. A few of the results of this 
paper have already been presented in a short paper [26]; the aim of the present paper is to show in detail how the 
RSRG method applies to such disordered spin models and to explore more of its consequences. Although we will give 
here a detailed discussion of the RSRG for the spin model, we will rely on Ref. [28] for many results about the single 
particle diffusion aspects of the problem; these we will only sketch, referring the reader to [28] for details. 

As for the single particle problem, the RSRG method allows us to compute a great variety of physical quantities, 
remarkably, including even some which are not known for the corresponding pure model (e.g. the domain persistence 
exponents 5 and ip). This provides another impetus for the study of the random models. Using the RSRG we also 
obtain the equilibrium behavior which corresponds to a well defined scale at which the decimation is stopped. 

The RSRG method is closely related to that used to study disordered quantum spin chains [30-35]. The crucial 
feature of the RG is coarse-graining the energy landscape in a way that preserves the long time dynamics. Despite its 
approximate character, the RSRG yields for many quantities asymptotically exact results. As in [32] it works because 
the width of the distribution of barrier heights grows without bound on long length scales, consistent with rigorous 
results of [27] . It is interesting to note that an exact RG has also been applied to the problem of coarsening of the 
pure ID soft-spin Ising $ 4 model at zero temperature for which persistence exponents have been computed [36,37]. 
Extensions to higher dimensions of the RSRG method for equilbrium quantum models have recently been obtained 
[38] ; these introduce hope that other dynamical models could be studied in higher dimensions as well. 

The outline of this paper is as follows. In Section II we define the RFIM and spin glass models and their dynamics 
and in Section III we summarize some of the main results and the physics involved. In Section IV we explain how the 
RSRG method is applied to the RFIM, give the corresponding equations and fixed point solutions and discuss the 
properties of the asymptotic large time state. In Section V we derive exact large time results for single time quantities 
such as the energy, the magnetization and the domain size distribution. In Section VI we obtain the time dependent 
single time two spin correlation functions both with and without an applied field. In Section VII we compute the two 
time spin autocorrelation as well as the two-point two-time correlations. In Section VIII we study aging phenomena 
which necessitate considering rare events which dominate all moments of the thermal (truncated) correlations, which 
we obtain, as well as the response to a small uniform magnetic field applied after a waiting time t w as in a typical aging 
experiment. We also compute the fluctuation dissipation ratio X(t,t w ). In Section IX we compute the persistence 
properties — the probabilities of no changes — of a single spin, of the local magnetization and of domains. Finally 
in Section X we obtain the distribution of the total magnetization and of the number of domains for large finite 
size systems. A brief discussion of the possibilities for experimental tests of the predictions and speculation on the 
applicability of the some of the general features to higher dimensional systems are presented at the end of the paper. 

Various technical results are relegated to appendices. In Appendix A the convergence towards the asymptotic state 
is shown. In Appendix B the time dependent single-time two-spin correlation functions in an applied field is derived 
while in Appendix C, the two-point two-time correlations are computed in detail and in Appendix D some of the 
finite size properties are analyzed. 



II. MODELS AND NOTATION 



A. Random Field Ising Model 



1. Statics 



We consider the random field Ising chain consisting of N spins {Si — ±l}»=i,jv with Hamiltonian : 



n=N-l 



i=N 




(1) 



with independent random fields with identical distribution whose important moments are: 



H=h t 



(2) 
(3) 
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so that H can be considered as a uniform applied field and g is the mean-square disorder strength; here and henceforth 
we denote averages over the quenched randomness — usually equivalent to averages over different parts of the system 
- by overbars. 

The statics and the dynamics of the random field fsing model (RFIM) can be studied in terms of domain walls 
living on the dual lattice: on the bonds + 1)} will be indexed by their left-hand site i. A configuration of spins 
{Si} can be represented as a series of "particles" A corresponding to domain walls of type (+|— ) at positions a\, ..(In a , 
and "particles" B corresponding to domain walls of type (— 1+) at positions b\, ..&/v B ; these must of course occur in 
an alternating sequence ABABAB.... [The relation between the numbers Na and Nb of domain walls depends on 
the boundary conditions : for instance in a system with periodic boundary conditions, Na — Nb, while in a system 
with free boundary conditions, \Na — Nb\ < 1-] 

It is very useful to introduce the potential felt by the domain walls: 

X 

V{x) = -2Y,hi (4) 

i=l 

and rewrite the energy of a configuration {Si} in terms of the positions of the set of domain walls as: 

N A N B 

H = H ref + 2J(N A + N B ) + J2 V M - E V ^ ( 5 ) 

a—l a—1 

where H re f is the energy of the reference configuration where all spins are (+) : 

n—N 

H re f = -J(N-l)- hi ( 6 ) 

i=l 

Each domain wall costs energy 2 J; the domain walls A feel the random potential +V(x); whereas the B walls feel the 
potential — V(x), i.e. the two type of domain walls feel opposite random potentials. 

Let us first recall some known features of the statics. In the absence of an applied field (H = 0) the system is 
disordered at T — and contains domains with typical size given by the Imry Ma length: 



4J 2 

-<IM ~ 

9 



(7) 



obtained from the following simple argument: the creation of a single pair of domain walls (A, B) a distance L apart 
costs exchange energy 4J independent of L. But the random potential has typical variations \V(x) — V(y)\t yp ~ 
\/kg\x — y\, and thus the typical energy that the system can gain on a length L by using a favorable configuration 
of the random potential is of order \/4gL. The two energies become comparable for L ~ Ljm and for L > L/m, 
it becomes favourable for the system to create domain walls. Thus the ground state will contain domains of typical 
size Ljm- One should note the difference between the case of e.g. bimodal distributions (for which the ground state 
can be degenerate) and continuous distributions for which it is non-degenerate; we will generally consider continuous 
distributions for simplicity. 

At positive temperature without a random field, the thermal correlation length Lt ~ exp(2J/T) gives the typical 
size of domains in the system in equilibrium at temperature T. But in the presence of the random field the equilibrium 
state at finite temperature is dominated by the random fields and still given by the Imry Ma picture provided: 

1 < Ljm < L T (8) 

which is the regime studied in the present paper. 

In presence of a uniform applied field H > the system at low temperature will contain domains of both orientations 
but with different typical sizes leading to a finite magnetization per spin m(g,H) < 1. 



2. N on- equilibrium Dynamics 

In this paper we will study the magnetization non-conserving dynamics of the RFIM model (1) starting from a 
random initial condition at time t = corresponding to a quench from a high temperature state. 
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Although our results are independent of the details of the dynamics (provided it satisfies detailed balance and is 
non-conserving) let us for definitencss consider Glauber dynamics where the transition rate from a configuration C to 
the configuration Cj obtained from C by a flip of the spin j, is: 



-0AE 



e /3AE _|_ e -f3AE 



(9) 



which satisfies the detailed balance condition and where AE = 2JSj(Sj-i + Sj + \) + 2hjSj is the energy difference 
between the two configurations, which takes the following possible values in terms of domain walls: 



AE{ creation of two domain walls} = 4 J ± 2hj 

AE{ diffusion of one domain wall} = ±2hj 

AE{ annihilation of two domain walls} = —4 J ± 2hj 



In this paper we focus on the regime: 



{hi} < J 



(io) 
(ii) 

(12) 



(13) 
(14) 



in which (8) holds. The randomness will dominate on length scales that are sufficiently large so that the cumulative 
effect of the random field energies is greater than T; i.e. for scales 



I > h 2 /T 2 . 



(15) 



So that there will be only one basic length and time scale, it is simplest to consider {hi} <~ T. 

There will be a wide range of time scales during which the domains grow from initial sizes of O(l) to sizes of order 
Lim by which time the system will be close to equilibrium. Since the Imry-Ma length Liu is very large, we expect 
universality of the long time dynamics, as we indeed find. 



B. Ising spin glass in a magnetic field 

In this paper we also consider the ID Ising spin glass in a uniform field h: 

i=N-l N 

H = - Ji(Ji<Ji+\ - ^2 ha l (16) 

i=i i=i 

As is well known, in the case of a bimodal (± J) distribution with equal probabilities for either sign, ID Ising spin 
glasses in a field are equivalent via a gauge transformation to random field ferromagnets. More precisely, setting 
Ji = Jei = ±1 for i = 1, ..N — 1, and defining a\ = Si and Oi = e\..ei-\Si for i = 2, ..N, the gauge transformation 
gives the new Hamiltonian 

i=JV-l N 

H =- E JSiSi +1 -Y,hiSi (17) 

i=l i=l 

where hi = he\..€i-\. Since the q are independent random variables taking the values ±1 with probability 1/2, the 
fields hi are also independent random variables taking the values ±/i with probability 1/2. Hence the new Hamiltonian 
describes a (±/i) random field Ising model with H = and g = h 2 in (1). 

The physical interpretation is as follows. The spin glass chain in zero field (h = 0) has two ground states, given by 
±cr| ' , where af^ = e\..ei-\, which correspond, via the gauge tranformation, to the pure Ising ground states sf^ = +1 
and S\ 0) = -1. In the presence of a field h > 0, the ground state of the spin glass is made out of domains of cither zero 
field ground states. Their typical size is thus given by the Imry-Ma length Ljm — ^r- These domains correspond 
to the intervals between frustrated bonds since at the position of a domain wall one has JiGi<Ji+\ = JSiSi+i < 0. 
Similarly each domain has magnetization: 

I E a i 0) \ = l\ E hi\ = ^\V(a a )-V(b a )\ (18) 

?G domain z£ domain 
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which is thus proportional to the absolute value of the corresponding barrier of the Sinai random potential, a property 
which will be used below. 

Via the gauge transformation, the non equilibrium dynamics (e.g. Glauber dynamics) of the spin glass in a field 
starting either from (i) random initial conditions (corresponding to a quench from high temperature) or from (ii) the 
pure ferromagnetic state (obtained by applying a large magnetic field that is quickly reduced to be h << J) corresponds 
to the non equilibrium dynamics (e.g. Glauber) of the RFIM also starting from random initial conditions. Thus in 
the regime h ~ T <C J many of the universal results obtained for the RFIM will hold directly for the spin glass in a 
field. Thus we will study the two models in parallel in the present paper. 

An important complication in applying the results to ID spin glasses is that the strengths of the exchange couplings 
will in general be random. This provides an extra random potential for the domain walls which has, however, only 
short range correlations. The conditions for there to be a wide regime of validity of the universal coarsening behavior 
found here is that the distribution of exchanges be narrow: 

\J\max |^|mm < \ J\ ■nun • 

(19) 

ideally, with |J| 

max \J\min ^ h ot smaller. Note that the equilbrium correlation length will be set by \ J\min for 
small h as the domain walls can easily find positions with | J| = |J| m j n that are near to extrema of the potential 
caused by the random fields. 

Having defined the models, before turning to the details of the calculations, we now briefly decribe the important 
physics and summarize our main results. 



III. QUALITATIVE BEHAVIOR AND SUMMARY OF RESULTS 

As for any Ising system, the static and dynamic properties of random spin chains can be fully described in terms 
of domain walls between "up" and "down" domains. The crucial simplification in ID is that the domain walls, which 
come in two types A = (+|— ) and B = (— 1+), are point objects. As discussed above, a random field induces a 
potential that the walls feel which has the statistics of a random walk so that its variations on a length scale L are 
of order \[L. The A walls tend to minima of the potential while the B walls tend to maxima. If they meet they 
annihilate but, on the time scales of interest here, the probability that a pair is spontaneously created, e~ 4J ' T can be 
ignored. 

As shown by Sinai [27], the motion of a single domain wall in such a random potential is completely dominated by 
the barriers that have to be surmounted to find low energy extrema. Since the time to surmount a barrier of height 
b is of order e b / T ', and to find a minima a distance L away a barrier of order b <~ h\[L will have to be overcome, the 
typical distance a wall moves in time t is only of order 1(f) ~ In 2 t. Although this motion is controlled by rare thermal 
fluctuations that take the wall over a large barrier, the position of a wall that started at a known point can, at long 
times, be predicted with surprising accuracy. This is because the width of the distribution of barriers to go distances 
of order L is as broad as the magnitude of the lowest barrier which enabled the wall to move that distance. Thus at 
long times, when the barriers become very high, the probability of going first over other than the lowest surmountable 
barrier that delimits the region in which the wall is currently, is extremely small. The position of a single wall at long 
times is thus determined by its initial position and by the height of the maximum barrier which it could surmount up 
to that time; this quantity, which we denote 

r = Tlnt (20) 

yields a well defined region that the wall can explore up to time t. The boundaries of the appropriate region that 
encompasses the wall's initial position are determined by T(t). The wall will be in local equilbrium in this "valley" 
on time scales of order t and thus tend to spend most of its time near the bottom of the valley. This behavior, as 
proved by Golosov [40], implies that the long time dynamics of a single wall is pseudodctcrministic: rescaled by the 
typical distance it has gone, In 2 t, the wall's position is asymptotically deterministic at long times. 

In a random field Ising chain, the two types of domain walls move in random potentials which are identical except 
for their sign. When walls meet, they annihilate. Not surprisingly, since each wall can move a distance of order In 2 t 
in time t and they cannot pass through each other or occupy the same position, the density of domain walls remaining 
at time t is simply of order 1/ In 2 t so that the correlation length starting from random (or short-range correlated) 
initial conditions grows as 



G 



£(t)~r 2 /5~ln 2 t. 



(21) 



This is in sharp contrast to the much faster power-law growth of the correlation length in most non-random systems, 
for example, the £(t) ~ \ft for Ising systems with non-conserved order parameter. 

The time-dependent correlation length can be more precisely defined from the average non-equilibrium correlation 
function. This, and most other non-equilibrium properties, are scaling functions of the ratio of lengths to powers of 
log-times, as is characteristic of "activated dynamic scaling" [41] which occurs in many random systems. We find that 
in zero applied field the average equal-time correlation function behaves as 

Jmm . = - (22) 



whose Fourier transform is at large time: 

+00 



tanh 2 (7iV^(t)/2) 



(23) 



with 3? denoting the real part. These results obtain until the non-equilibrium correlation length 



reaches the equilibrium correlation length 



n T 2 ln 2 t 

«*> = ^ (24) 



8</ 2 2 

£,eq = ~ o~ = —oLlM (25) 



It is intriguing that the form of the non-equilibrium correlations in the universal scaling limit are identical to those 
of the equilibrium correlations, the only difference being the correlation length. 

A remarkable property of the coarsening in the RFIM chain was conjectured in [26]: the positions of the set of 
domain walls and hence all of the correlations in the non-equilibrium state are asymptotically deterministic at long 
times. This means that while the evolution of the domain wall structure is only logarithmic in time, the domain wall 
configurations of two runs following quenches using the same sample converge to each other much more rapidly: as a 
power of time! More precisely, the probability that the spin configurations of the two runs measured at the same time 
after the quench differ substantially in a region of size of order £(i) decays to zero as a non- universal power of time, 
becoming negligible well before the system is anywhere near equilibrium. [Strictly speaking, because in rare region, 
as discussed later, fluctuations of the position of an individual or neighboring pair of domain walls will exist, one will 
need to do a certain amount of time averaging for this pseudodeterminism to become most evident] 

The asymptotic determinism will occur even if the initial conditions were macroscopically distinct in the two runs: 
for example, if one was quenched from a high temperature state in zero applied field, and the other from a high 
temperature state with a small net magnetization caused by a uniform applied field. This pseudodeterminism in a 
system with many degrees of freedom is a dramatic effect; it implies that history dependence can, even under strongly 
non-equilibrium conditions, be weak; one might thereby be fooled into thinking that such history independence implies 
equilibrium. 

While the equal-time correlations during coarsening converge to an equilibrium-like form, two-time quantities show 
interesting history dependence, generally depending on both times rather than just the time difference as they would 
do in equilibrium. There are typically three regimes: the later time, t, much longer than the earlier, t', the two 
times of the same order, and the difference between the two times, t — t' much smaller than either time. The scaling 
variables in these three regimes are, respectively, hit'/ hit, t'/t, and ln(t — t')/lnt. In the first and third of these, 
the condition for asymptopia, that the scaling variable is small, is very difficult to attain in practice, thus knowing 
the full form of the scaling functions is essential for analyzing experimental or numerical data. Note that only in the 
last of these regimes should one expect to find equilibrium-like behavior characteristic of the local equilibrium that is 
being probed. 

In this paper we compute a variety of correlations and response functions that illustrate some of the interesting 
non-equilibrium behavior. The simplest of these is the temporal autocorrelations of a single spin. With no applied 
field, the average autocorrelation function is found to decay as a power of the ratio of the correlation lengths at the 
two times: 
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(S x (t)S x (t>)) ~ (|^) (26) 

for t > t' with the exponent A = 1/2. In a three dimensional spin glass, the autocorrelation function determines the 
non-equilibrium decay of the magnetization after a large uniform field is turned off below the transition temperature 
[5]. There are complications here due to the fact that a small applied field needs to be left on to provide the 
randomness, but the basic physics is similar. 

The single spin autocorrelation function provides some information on how much memory the system has of its 
earlier history. More information is provided by the two-spin two-time correlation function (S x (t)S y (t))(S x (t')S y (t')) 
which converges to a scaling function of (x — y)/t;{t) and £(t')/£,(t) whose Fourier transform in (x — y), D(q,t,t'), we 
compute exactly, to our knowledge the first such computation for a physical model pure or random. For an ideal spin 
glass in which there are no correlations between the positions of positive and negative exchanges and the distribution 
of the exchange strengths is symmetric in J — > — J, this correlation function, D, is the average over the randomness — 
or equivalcntly a positional average over the regions being probed by the scattering — of the product of the magnetic 
scattering intensity at time t and wavevector k and that at time t' and wavevector q — k. 

Other properties associated with "aging" can also be computed exactly; we study the thermal fluctuations around 
the configurations at two different times and the dynamic linear response to a uniform magnetic field that is turned 
on after waiting for some time for the system to equilibrate. In equilibrium, these are related by the fluctuation 
dissipation theorem but we find, as expected, that this relation generally fails except in the limit that the time 
difference is much less than the waiting time. The behavior we find is, however, rather different from that found in 
mean-field models [45] and we discuss the contrasts between these results. 

In random systems controlled by zero temperature fixed points, such as is the case in the regime of scales studied 
here, thermal fluctuations and linear response functions are both dominated by rare spatially isolated regions of the 
system — although which regions dominate depends on the time scales and properties of interest. The study of these 
thus involves the corrections to the deterministic approximation to the dynamics that led, as discussed above, to 
exact asymptotic results for many other quantities. The dominant events are rare by a factor of, typically, 1/lnt, but 
nevertheless still lead to universal results. 

"Persistence" properties provide another probe of how much memory a system retains of its initial configuration, 
for example, what the probability is that a spin has never flipped. Surprisingly, this and other related quantities — 
including some which are not known in pure systems — can also be computed exactly for the RFIM chains. 



IV. REAL-SPACE RENORMALISATION METHOD 

The coarsening process taking place in non equilibrium dynamics of the RFIM starting from random initial condi- 
tions can be thought of as a reaction-diffusion process in a ID random environment (Sinai landscape) for the domain 
walls. Indeed, in the regime T <C J considered in this paper, the creation of pairs of domain walls is highly suppressed. 
The dynamics of the domain walls is thus dominated by the random field as follows. The domains walls A quickly 
fall to the local minima of the random potential V(i) = Vi, whereas the B walls quickly move to the local maxima of 
Vi. Then they slowly diffuse by going over barriers in opposite Sinai potentials ±Vi. When an A and a B meet, they 
annihilate preserving the alternating ABAB... sequence. 

We can thus use the real-space renormalisation procedure introduced in [28] to study the Sinai diffusion of a single 
particle, and extend it to take into account the annihilation processes. In the single particle case this procedure was 
shown to be asymptotically exact at long length and time scales. Here it is also expected to be asymptotically exact at 
large time, since as in the Sinai case, we find that the the effective distribution of the random barrier heights become 
infinitely broad in the limit of large scales. Thus the motion over the barriers becomes more and more deterministic 
at long times. 

A. Definition of the real space renormalization procedure 

We briefly outline the renormalization procedure, detailed in [26,28], for the diffusion of a particle in the Sinai 
landscape V,. Grouping segments with the same sign of the random field, one can start with no loss of generality from 
a "zigzag" potential Vi where each segment ("bond") is characterized by an energy barrier Fi = |Vj — and a 
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(b) 




li h 1 3 

FIG. 1. (a) Energy landscape in Sinai model (b) decimation method: the bond with the smallest barrier F m i n = F2 is 
eliminated resulting in three bonds being grouped into one. 

length li. From the independence of the random fields on each site, the pairs of bond variables (F, I) are independent 
from bond to bond and are chosen from a distribution P(F, I) normalized to unity. 

The RG procedure which captures the long time behaviour in a given energy landscape, is illustrated in Fig. 1. 
It consists of the iterative decimation of the bond with the smallest barrier and hence the shortest time scale for 
domain walls to overcome it [26] [28]. This smallest barrier, say, F2, together with its neighbors, the two bonds 1,3, 
are replaced by a single renormalized bond with barrier 

F' = Fi - F 2 + F 3 (27) 

and length 

l' = h+h + l 3 . (28) 

This defines a renormalized landscape at scale T where all barriers smaller than T have been eliminated. 

Since the distribution of barriers is found to become broader and broader [32] an Arrhenius argument implies that 
the diffusion of a particle becomes better and better approximated at large time by the following "effective dynamics" 
[26,28]. The position x(t) of a particle that started at xo at t — coincides with - or is at least very close to - the 
bottom of the renormalized bond at scale 

T = Tlnt (29) 

which contains xq. Note that we choose time units so as to set the microscopic (non- universal) inverse attempt 
frequency to unity [26,28]. This RG procedure is thus essentially decimation in time. Processes that are faster than 
a given time scale are decimated away and assumed to be in local equilibrium. 

In the presence of domain walls of type A and B, we must keep track of both the diffusion and the possible reactions 
of the domain walls that occur during the decimation. Upon the decimation of bond 2 (see figure 1), there are four 
possible cases, illustrated on figure 2, according to whether or not there is an A at the bottom of bond (2) and whether 
or not there is a B at the top of bond (2). If (i) there are no A, no B, then one simply renormalise the bond and 
nothing happens to the domain walls. If (ii) there is an A but no B then A goes to the bottom of the new renormalized 
bond. If (iii) there is a B and no A then B goes to the top of the new renormalized bond. And if (iv) there is both an 
A and a B then the two domain walls meet and annihilate upon decimation. This annihilation will occur at a time of 
order e r / T determined by the barrier of the decimated bond: by assumption, T at this scale. 

The only truly new process compared to the single particle diffusion is the case (iv). The above RG rule is consistent 
with the large time dynamics in this case for the following reason. Since A and B domains diffuse with independent 
thermal noises, the time it takes for them to meet is again Tint — F (accurate on a log scale) . Indeed the probability 
that they meet at a point at potential V is ~ exp(— (F — V)/T) exp(— V/T) = cxp(— F/T) (as can be seen from 
considering the equivalent 2d diffusion problem (of the pair of wall positions) with an absorbing wall at x + y = I). 

The second difference from the procedure in the case of the single particle, is that it must be stopped when the 
renormalisation scale T reaches the energy cost of creation of a pair of domain walls: 
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r = Tj = 4j 



(30) 



Indeed beyond this scale Tint > Tj one must take into account creation of pairs of domain walls. As will be shown 
below, the state at T = Tj gives the final equilibrium state. 




FIG. 2. The four possible cases of the evolution under decimation of the middle bond of the domain walls, given the constraint 
of alternating sequence of domains. A and B denote respectively the domain walls (+| — ) and (— 1+), denotes no domain wall 
present at the top or bottom of the renormalized bond, while C represents either a B or and a D represents either a A or 0. 



B. RG equations and statistical fixed point of the landscape 

Independently of whether the bonds contain domain walls (A or B) or not, which is studied in the next section, 
one can study the evolution under RG of the landscape. Since the RG rules preserve the statistical independence 
of the variables (F, I) from bonds to bonds, it is possible to write closed RG equations for the landscape, i.e. for 
Pp(C = F ~ T, l) and Pf (£ = F — V, I) which denote the probabilities that a ± renormalized bond at scale T has a 
barrier 

F = r + c>r (3i) 

and a length I, each normalized by J °° d£ J °° dlP^(C, I) = 1. 

(dr - d ( )P±({,l) = P r T (0, .) P r ± (., .) * u P r ± (., .) (32) 

+Pr(C,0 [°° dl'(P±(0,l')-P*(0,l')) (33) 
Jo 

where denotes a convolution with respect to I only and with respect to both £ and I with the variables to be 
convoluted denoted by dots. As discussed in [32] [28], the solutions of these RG equations depend on an assymetry 
parameter S defined as the non vanishing root of the equation 



e- iSh = 1 (34) 

which reduces in the limit of weak bias H to 

8 - ?~ ( 35 ) 
with 5 = in the absence of a uniform applied field. Our results arc valid for long times as long as 
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ST < 1 (36) 

For large T, the Laplace transform of the distributions take the following form, in the scaling regime of small S 
and small p with 5T fixed and pT 2 fixed [32] 



f 

Jo 



dlP±(U)e~ ql = C/±(^-)e- c "r(*) (37) 



(p) = v / p + <5 2 coth pVp + <5 2 ] t 5 



r ; sinh [rvp+^] 



This exact knowledge of the renormalized landscape will be used below to extract physical quantities for the spin 
models. 



C. Convergence towards "full" states: asymptotic determinism 

Even armed with the statistical properties of the renormalized landscape, we still have to determine the long time 
distribution of the occupation of the extrema of this landscape by A and B domain walls. At first sight this seems a 
difficult problem. Indeed the positions of A's and B's can be correlated over many bonds of the renormalized landscape 
since there are a priori empty maxima and minima and the domain wall positions must respect the alternating 
constraint ABABAB. 




FIG. 3. Full state in the renormalized landscape; as we show, this corresponds to the state at large time (see text). Each 
top and bottom of bonds are occupied by respectively a B and A domain wall. The correponding spin orientations are also 
indicated for the RFIM. Decimation under increase of F corresponds to annihilation of the domain with the smallest barrier. 
For the spin glass the position of the domain walls corresponds to the frustrated bonds. 

However, the RG analysis becomes simple if the system reaches at some stage a "full" state which has one A wall 
at each minimum and one B wall at each maximum of the renormalized landscape as illustrated in figure 3.. It is easy 
to see that such a "full" state is preserved by the RG procedure. Also note that this "full" state would obtain from 
the beginning if, for instance, the initial condition were complctly antifcrromagnetic. 

Generally, we consider random initial conditions and thus the initial state is not a "full" state. However, one can 
show that the system converges exponentially in T towards a "full" state, as we now discuss. 

The renormalisation procedure for the coarsening process of the RFIM has the following important property : the 
configuration for the spins at scale T depends only on the renormalized landscape at T and on the initial configuration 
of the spins (e. g. equilibrium at high temperature before the quench at t = 0). In particular it does not depend on 
the initial landscape or — except occasionally at early times before the barrier distribution becomes very broad — 
on the whole history of the reaction-diffusion processes of domain walls. 

Let us first consider one ascending bond with extremities (x, y) of the renormalized landscape and assume that there 
were initially n domain walls in the interval (x, y). Neglecting for the time being the influence of the two neighboring 
bonds, it is easy to see that the state at T is determined only by n and does not depend on the order in which the 
reactions between domain walls have occured: indeed, it is determined by the parity of n and by the nature of the 
domain wall closest to the bottom end. There are several cases 

i) If n = 0, then the final state of the occupation of the endpoints of the bond is (0, 0). 

ii) If n is odd, and if the domain wall closest to the bottom is of type A, then the final state is (A, 0) 

iii) If n is odd, and if the domain wall closest to the bottom is of type B, then the final state is (0, B) 

iv) If n is even with n > 2 and the domain wall closest to the bottom is of type A, then the final state is (A, B) 
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v) If n is even with n > 2 and the domain wall closest to the bottom wall is of type B, then the final state is (0, 0). 

Of course, to obtain the real occupation of the top and the bottom ends of one renormalized bond, one needs 
to consider what happens also on the two neighboring renormalized bonds and to compute the probabilities of the 
various states, taking properly care of the alternating constraint ABAB. This is done in Appendix A. It is found that 
the crucial feature is that in order for a bond not to have both ends occupied by domain walls, either it or one of 
its neighboring bonds must have had no domain walls on it initially. Since the bonds tend to become progressively 
longer with time, the chance of this occuring drops rapidly with increasing T. This will be true even if the positions of 
the domain walls have some local correlations, in particular for the case of an initial state that corresponds to a high 
temperature configuration in a small magnetic field so that there is an initial magnetization and the typical distances 
between a neighboring A and B wall will depend on which of the two is on the left. 

We thus find that the system converges towards the "full" state of the renormalized landscape exponentially fast in T, 
with a non-universal coefficient that depends on the initial concentration of domain walls, on the inital magnetization, 
and on the strength of the randomness. This corresponds to a power law decay in time, as t~ n with n non-universal, of 
the probability that a maximum or minimum of the renormalized landscape at time t is unoccupied. The positions of 
the full set of domain walls at long times are thus asymptotically deterministic and independent of the initial conditions 
provided these have only short range correlations. Concretely, this asymptotic determinism can be characterized by 
the typical mean distance between "missing" domain walls, i.e., deviations from the deterministic full state. The 
distance between these "errors" will grow exponentially in Y and hence as a power of time. The only exceptions 
to this are associated with fluctuations in the (nominally) full states: because of the rare times or configurations in 
which domain walls are either about to annihilate or are spending substantial fractions of their time in more than one 
valley; these fluctuation induced "errors" we will study in detail later; they decay far more slowly with time than any 
persistent initial condition induced differences between runs. 

Since all the universal quantities that we will study have length scales that grow much more slowly - typically 
logarithmically — in time, or probabilites of occuring which decay much more slowly with time - typically as powers 
of 1/r ~ 1/ In 2 t - we restrict consideration in what follows to the analysis of full states. Note that our results for the 
biased case in which there is a small uniform field applied in addition to the random field do yield power law growth 
of the correlation length with time, but in the regime of validity of these results, the power law is very small and 
again the effects of "missing" domain walls are negligible at long times. 

D. Convergence towards equilibrium 

As mentioned earlier, the RG procedure must be stopped at 

r = Tj = AJ (38) 

since at this scale one must start to take into account the energetic benfits of creation of pairs of domain walls. At 
scale r = Tj the typical domain size is the Imry-Ma length Ljm — 4^J 2 /g, and the energy cannot be lowered further 
by any process in the full state. Indeed, moving the domains walls without changing their number cannot lower the 
total energy since domain walls already occupy all tops and bottoms of the renormalized potential. Decreasing the 
number of domain walls by two can also not lower the energy since the gain is 4 J while the loss due to the random 
field for walls separated by a bond of barrier F is F > T = AJ. Similarly, adding two walls the cost is AJ and the 
gain is F < T = AJ since the only positions they can occupy are by definition separated by a barrier < T which 
has already been decimated. Thus if the renormalization is stopped at T,/, in the small field, low T scaling limit 
the configuration of the walls corresponds precisely to the ground state and, up to negligible thermal fluctuations, 
to the thermal equilibrium state. Thus we are able to compute equilibrium properties straightforwardly from the 
renormalization group analysis. 

Thus the RG allows one to study the approach to equilibrium starting from any initial condition characterized by 
typical domain sizes Loe 2J / T ° <C Ljm, where Tq represents the temperature before the quench. As explained above, 
under these conditions the relaxation towards equilibrium always takes place by diffusion and annihilation of domain 
walls before any domain creation can occur. 

E. Approach to equilibrium from more ordered initial conditions 

A qualitatively different type of evolution towards equilibrium takes place in the other limit (not studied here): 
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A J 2 

U ~ e 2J/To » L IM = ±L (39) 

still imposing h <C To . This is a regime of initial temperatures low enough that the initial density of domains is very 
small. In that case a very different relaxation process towards the same equilibrium state discussed above takes place: 
for r < Tj = 4J well separated domain walls diffuse independently with very rare annihilations. When T reaches Tj 
many domain walls are "suddenly" - within a factor of two or so in time - created and the large initial domains break 
into many smaller ones of size Ljm- Thus, the relaxation is more abrupt in this case than in the more interesting one 
studied in the present paper. 

V. ENERGY, MAGNETIZATION AND DOMAIN SIZE DISTRIBUTION 

From the knowledge of the fixed point for renormalized landscape (37) and the fact that the system reaches the full 
state (exponentially fast in T) where each top is occupied by a B domain and each bottom by a A domain, we can 
immediately compute several simple quantities. We will compare these results with the existing exact results known 
in the statics. 

A. RFIM without applied field 

Specializing (37) to the case S = in the absence of an applied field the fixed-point of the RG equation is [31,32] 

P*(„, A) = LT;\ (-^L= e-Vicothv^ {m) 

where we have introduced the dimcnsionless rcscalcd variables for barriers 

r ] =(F-T)/T (41) 

and for bond lengths 

A = 2g^ (42) 



1. Number of domain walls per unit length 

Thus for large T = Tint, the average bond length Ir, equal to the average distance between two domain walls, 
behaves as 

7 r = — r 2 = — T 2 \n 2 t (43) 



4ff 4 3 

The number of domain walls per unit length decays as 



"(*) ^ ( 44 ) 
T 2 In 2 1 



up to time t eq ~ exp(4J/T) at which equilibrium is reached and: 

1 



n(t eq ) = n eq = (45) 
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2. Energy density 



The energy per spin as a function of time (i.e of T = Tint) is simply given by 



E r ~-J + n r (2J-^<F> r ) (46) 

=- J+ w^ {2J - Tlnt) (47) 

where < F >r denotes the averaged barrier at scale T. This formula holds up to time t eq <~ exp(4J/T) where the 
ground state energy at equilibrium is reached: 

E as--J-{j (48) 

Since we consider the regime g -C J, this result is expected to be exact to first order in g/J. It does indeed agrees 
with the exact result of [20] concerning the bimodal distribution, expanded to first order in g. To obtain higher orders 
in the expansion in g/J one would need to compute within the RG higher orders in a 1/T expansion. 

Note that the entropy per spin at T = computed in [20] for ±h distributions originates from degenerate configu- 
rations occuring from short scales and is thus non universal. If the distribution is continuous we expect S ~ T from 
short scales, also non universal. 



3. Distribution of lengths of domains 

Since the bond lengths in the renormalized landscape are uncorrelated we obtain the result that the lengths of the 
domains in the RFIM (both in the long time dynamics and at equilibrium) are independent random variables. [Note 
that this is different from the exact result for the dynamics of the pure Ising chain obtained in [39] .] 

Moreover, during the coarsening process, the probability distribution of the rescaled length A = ^ = jSh? t * s 
obtained as: 



CO s v CO ^ 

E ( n+ ^<- l ^ x{n+hr = 7^ E (-D"v + >-^ 

i — — rvi ^ ' " 'm — — ("v— i 



(50) 



The distribution of the length of domains at equilibrium is also given by P(X eq ) where: 

2gl 



^eq - - o7 ( 51 ) 

1 7 ZLr 



B. Spin glass in a field 

Using the gauge transformation described in Section (II B) the above results readily apply also to the spin glass in 
a field. Let us recall that a domain in the RFIM corresponds in the SG to an interval between two frustrated bonds. 
Using the above expressions for the zero applied field RFIM and replacing g — > h 2 , we thus obtain the averaged size 
of these domains from (43), their number per unit length from (44) and (45), and their distribution of lengths from 
(49) and (51). 

The distribution Pr(M) of the magnetization M — \ J2iedomain a i\ °f eacn domain is obtained from the distribution 
of barriers as: 

1 M - Mr 

Pr{M) = Mr CXP( M7^ (M ~ Mr) (52) 

M T = L (53 ) 
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with r = Tint. In equilibrium the same result holds with Mr — > 2 J /h. Note that since this variable is proportional to 
the barrier, there are no domains of magnetization smaller than Mr. Similarly the joint distribution of magnetization 
and length is given by (40) in Laplace transform and rescaled variables 77 = M ~^ T and A = 2h 2 l/T 2 , with again the 
property of statistical independence of domains. 

Finally, the energy per spin is given by the expression (47) replacing g — > h 2 . 

Note that the results here are, strictly speaking, restricted to the case in which the mapping to a random field 
Ising model is exact: the case in which all of the exchange interactions have the same magnitude and only differ in 
sign. Nevertheless, in the more general case with a distribution of | J|'s, the universal aspects of the non-equilibrium 
behavior will be the same as long as there is a non-zero lower bound to this distribution, |J| m j„. However at times 
longer than Tint — \ J\ m in, domain walls will no longer necessarily be annihilated; whether they are or not will depend 
on the local J as well as on the renormalized potential. This will, of course, also affect the equilibrium positions of 
domain walls, but because there will always tend to be weak exchanges near to the extrema of the potential caused 
by the random fields, the changes in the positions of the walls will be negligible on the scale of the correlation length. 



C. RFIM in an applied field 

In a similar manner to the above analysis, we obtain results when a small uniform field H is applied. The solutions 
(37) of the RG equations now depend on the parameter 5 defined as the non vanishing root of (34), equal to 5 = H/ (2g) 
in the small field limit for which our results will be asymptotically exact. 



1. Number of domain walls and magnetisation 

The averaged sizes of domains (respectively with spins oriented against and along the field) are found both to 
grow with time as: 

7- = r 2 J-(i-^lo (54) 

-+ 9 1 , „ sinh(7) N , „ 

ll=Y 2 — {e l ^-1 55 

where 7 = TS = 5T\nt and T = Tint. In the long time limit it is: 

l ^^2H (56) 
It « ^t™' 9 (57) 

Note that the fact that the length of domains with spins in the opposite direction from the applied field still grow on 
average, is simply due to the fact that the smallest ones (with barriers smaller than Tint) keep being eliminated. 
The number of domain walls per unit length thus decays as 

2 S 2 

n(t) ~ — = 4g 5 (58) 

J++J- y sinh 2 («Tnt) V ' 



and the magnetisation per spin grows as 



m(t) = ^ — =r = M\—T\nt] (59) 
7r + 7 r 2 5 

A4[7] = coth(7)- , 7 (60) 
sinh(7) 2 

The function A4[j] starts as ^[7] ~ §7 for small 7 and goes exponentially to A4 = 1 for large 7. 

These results hold up to time t eq ~ exp(4J/T) at which equilibrium is reached. The number of domain walls per 
unit length in the equilibrium state is thus: 



15 



s2 



with averaged sizes: 



L 



The equilibrium magnetisation per spin m eq is: 



sinli (4c5J) 



2J 



7+ ~ _9_ p 4JH/g 



/ \ i r ^T 

m e(? = m(i e(? ) = M[— — J 



(61) 

(62) 
(63) 

(64) 



2. Energy per spin 



The energy per spin as a function of time (i.e of T = Tint) is simply : 



E r = - J — H + n r (2J - - < F >_ 



-J-H + 



A5 2 g 



sinh^r) 



2j -K r+ ^ (i 



-2T5 



(65) 
(66) 



where < F denotes the averaged barrier of bonds against the field at scale T. This obtains up to scale T eq = 4J, 
at which equilibrium is reached with ground state energy 



E, 



gs 



-J-H 



sinli (4J<5) 



(1 



-8.75 



(67) 



This result is compatible with the result Eq. (80) of Derrida-Hilhorst (DH) [19] obtained from studying products 
of random matrices, as can be checked with the correspondance arm = 2(5. In addition we obtain here the explicit 
scaling form in the small 8 limit with 5J fixed. We have checked that this scaling form is also consistent with the 
exact result (equation (8) in [20]) for the bimodal distribution at leading order in S. 



3. Distribution of domain lengths P^(l) 

As in the zero field case, the lengths of the domains are independent random variables. Their probability distribu- 
tions can be obtained from Laplace inversion of (37). They read respectively for ± domains: 

+oo 

p ± (o = E c »W s «W e "' s " (7) ( 68 ) 

where 7 = STlnt and the functions ^(7), s„(7) are given in equations (50,51,52,53) of [28]. 



VI. EQUAL TIME TWO SPIN CORRELATION FUNCTION (S (T)S X (T)} IN THE RFIM 

We can compute the disorder averaged two spin correlation function by noting that, in a given environment, the 
equal time two spin thermal correlation (So(i)S x (t)) equals +1 if the points and x are at scale V = Tint on 
renormalized bonds of the same orientation (i.e. both ascending or both descending), and equals —1 otherwise. 
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A. Zero applied field 



The average over the environments is, in zero field H = 0: 

oo 

(S Q (t)S x (t))=Y,(- 1 TQr ) (x) 



(69) 



n=0 



where (x) is defined as the probability that point x belongs to the bond n given that the point belongs to the 
bond of the rcnormalized landscape (see figure 4). We have the normalization X^^Lo Qr ( x ) = 1- For n = 0, the 
probability that x is on the same renormalized bond as is 



1 f 00 f la f l ° 

Q { °\x) = = dl P T (l ) d yi dy 2 5(x-(y 2 - yi )) 

tr JO JO Jyx 

and thus after rescaling the variables, we find 

Q { °\x)= qo (x = 2g^) 

where 

q (X) = 2 / d\ {Xo-X)P*(\ ) 
Jx 

in terms of the fixed point solution P*(A) given in (49) for the distribution of rescaled lengths A = 2gl/T 2 . 

l \ 1 2 ^ 



(70) 



(71) 



(72) 




FIG. 4. Renormalized landscape at scale F — Tint indicating the sites and x of two spins and the bonds between them. 
For n > 1, we have 



Q^\x) = q n (x = 2g 



with 



q n (X) = 2 f P*(A )P*(Ai)...P*(A n _i)P*(A„) 

•> 1/1,1/2, Ai>0 

(5 (X - ( yi + Ai + A 2 + ... + A„_! + y 2 )) 9(X - yi)6»(A„ - y 2 ) 



The Laplace transforms read 



(*) = / dXe- sX q (X) = - [P*(s) - 1 
Jo s v 



r°° i — p*(<t\ 

q n (s)= dXe~ sX q n {X) = 2{ ^) 2 [P* (s)]^ 1 for n>l 

Jo s 



and thus 



£(-i) n </ n (*) 



1 4 /1-P*(s) 



n=0 



s s 2 \l + P*(s) 



(73) 

(74) 
(75) 

(76) 
(77) 

(78) 



Using now the explicit solution P*(s) = 1/ cosh(y / s), we finally get 
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(So(t)S x (t)) = LT~^ x=2gi 



1 4 ,2/Vs 

-tami — 

s s 2 V 2 



(79) 



V 48 + 64(2n + l)V^ _ (2 „ +1)2 , 22g M 



n— — 00 



(2n+ 1)% 4 



with r = Tint. The leading long-distance behavior of the equal time correlation function is proportional to xe 
rather than a simple exponential. The Fourier transform of the spin-spin correlation function is simply 



+00 

(81) 



ta,nh 2 (iry / iq£{t)/2) 



with 3? denoting the real part. As explained previously, the renormalization procedure has to be stopped at T ~ Tj = 
4J, i.e t <~ f e g ~ exp(4J/T) where the equilibrium state has been reached. Thus equation (80) with r = Tj gives 
the mean equilibrium spin correlation function. The correlation length £ (t) is given by the decay of the (n = 0) term 
which dominates at large distances 

^ T 2 ln 2 t 

^ = m 

up to scale Tlnt eq = Tj = 4 J and we obtain the correlation length at equilibrium: 

8 «/ 2 2 

Zeq = €(teq) = = ~2 L IM (83) 

This formula is in agreement with the limit h <C J of the exact result for the equilibrium correlation of [21] in the 
case of a bimodal distribution (±h). 

Finally, note that the RFIM two-point correlation function (Si(t)Si +x (t)} also corresponds for the spin-glass to the 
following correlation function involving the zero-field T — ground-state cr- ^ 



(Si(t)S i+x (t)) = ^fM+l^K+xW) (84) 



B. Non-zero applied field 

In presence of an applied field H > 0, one defines C t ei,C2 (x) as the probability that sign[So(t)] = £\ and sign[S x {t)\ 
e 2 - One has: 



(So(t)S x (t)) = C++(x) + Ct~{x) - 2C+-(x) (85) 

We know already that 

7+ 

C++(x) + C+-(x) = Prob{sign[S (t)} = +1} = _ + r __ (86) 



Cr~(x) + C+-(x) = Prob{sign[S (t)} = -1} = (87) 



h 

7r +Ir 



and we can thus write the correlation function as 



It^ It 



(S (t)S x (t)) - (S (t)) (S x (t)) = 1 - 4C+~(x) - I (88) 

\l r + l r J 

Performing an analysis similar to the case of zero field, we obtain the Laplace transform 
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' d xe-«*Ct-(x) = (l-^( g ))(l-^r( g )) (89) 
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The Laplace inversion can be performed as in [28]. The correlation decays as a sum of exponentials, and the term 
which dominates the asymptotic decay gives a correlation length 

m = (90) 

2 gso (7) 

where r = Tint, 7 = ST In t and the function Sq (7) is defined by equations (50) and (52) in [28]. In particular the 
asymptotic behaviour for large 7 is: 

W) ~% ~ ■£ S t BT/a (91) 
Note that even for the equilibrium case (7 = STj — AS J), this correlation length: 

does not seem to have been obtained previously; it is very different from the correlation length of the truncated 
correlations - i.e. that of the thermal fluctuations - computed for bimodal distribution in [21] and discussed here in 
Section (VIII A). 



VII. AGING AND TWO TIME CORRELATIONS 



Some of the most interesting properties of random systems involve "aging" , the dependence of measured quantities 
on the history of the system, particularly on how long it has been equilibrated for. In this Section we study one of the 
fundamental properties which show the effects of aging: two time non-equilibrium correlations. As before, the system 
is quenched from a random initial condition at time t = and we study the aging dynamics at late times between 
time t' = t w - the waiting time - and t, both t' and t > t' being large. 

We first consider the autocorrelation function of a given spin, from which one can extract the autocorrelation 
exponent A. Note that since this is a single site quantity it applies directly to both the random field and spin glass 
problems. 



A. Spin autocorrelations in zero applied field 

We consider the autocorrelations of the random field Ising model in zero applied field, H = 0: 

C(t,i?) = (Si(t)Si(t>)) (93) 

Except at short times the system is in the "full" state and hence C(t,t') is simply the probability that the site i - 
which we take to be the origin - belongs at both time t' and time t to renormalizcd bonds with the same orientation. 
Thus, this quantity can in principle be obtained from the result in [28] for the probability P(x, t; x't'\0, 0) that a 
particle diffusing in a Sinai landscape starting at at time t = is at x' at t' and x at t. Indeed C(t, t') (since it is 
computed in the full state) is simply related to the probability that a particles postitions X{t) and X(t') have the 
same sign. However, it can also be obtained through a much simpler direct computation which we now present. 

We define Pjtjt(C) (resp. Pfr' (C) ) as the probability that the origin is on a descending bond at V , and is on a 
descending (resp. ascending) bond of strength £ at a later stage, T. The RG equations read 

(dr d c ) P±+ (C) = -2P*(0)P±+(Q + 2P r ^(0)P r ± (.) *< P±+ (.) + P±(.) * c P± (.)!?+ (0) (94) 

together with the initial condition 

W c) = U"*niP (95) 

P r -+,(C)=0 (96) 

We introduce the scaling variable: 
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Since for large T' , Pr> has reached its fixed point value (40), one has in terms of the rescaled variable rj = (/T: 

(ad a - (1 + V )d v + 1)P±+( V ) = 2 f dr ? 'e-( I '- f '')p±+(7 7 ') + V e^P^+(0) (98) 

Jo 

together with initial conditions at a = 1 

PS( , ) = i£i^TO = i (1+2 , )e _, m 

and f'aJiC'?) = 0. The solutions are 

P^ + (v)=l(Ai + +vBt + )e-" (100) 

with ^M^d? (101) 

Bt + =\±l (102) 

which obeys the normalisation condition J °° c£?7 (P+ + (n) + P~ + (n)) = |. Since with = we have Pa~(v) = 
we obtain 

/>oo 

C{t,t')= dr,(P++( V ) + p--(r,)-P+-(r,)-P-+(r,)) (103) 



= ^ + +^ + -^ + -B-+ = ±-^ (104) 
and thus the autocorrelation function of the RFIM in zero applied field at large times t > t' is 

4 /lnf \ 1 /lnt' x : 



<*<*»<'»- 5 UfJ- 5 J < 105 » 

In particular the asymptotic behavior for fixed t' is: 

imWT) oc (106) 

where Z(t) ~ In 2 i is the characteristic length of the coarsening in the RFIM and where the autocorrelation exponent 
is: 

A = \ (107) 
The auto-correlation being invariant under gauge transformations, we immediately obtain for the spin glass: 

fo(tMf)> = (Sim(t')) ~ with \ = \ (108) 

Note that this value of A saturates the lower bound of d/2 in contrast to the pure ID Ising case which saturates 
the upper bound of A = d [5]. 

B. Autocorrelations for the RFIM in an applied field 

In presence of an applied field H > the calculation is similar to the above ; it is detailed in Appendix B. Using 
the scaling variables: 
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1 = ST^^-T\nt (109) 

y = ST'^— Tint' (110) 
2g V ' 



for small H and the magnetization per spin. 

m(t) = A4(-y) = coth -y — 

sinh 7 



< Si(t) > =m(t) = M(j) = coth 7 - — ] — , (111) 



the result for the autocorrelation function is: 

< SiWSiit) >-< Si(t) > < Si(t') > = ((27 - i)M(i) + 7' coth 7 ' - 1) (112) 

sinh 7 



The long time asymptotic behaviour is: 



AHTlnt 



< Si(f)Si(t) >-< Si(t) > < Si(t>) > = g m(t>). (113) 

This result is valid, strictly, to leading order in ST. In general the exponent for the power law decay in time which 
occurs here and for other quantities in the presence of a uniform applied field will have 0(S 2 T 2 ) corrections. 

C. Two-point two-time correlation function for the spin glass 

In order to characterize the spatial aspects of aging dynamics in the spin glass we have computed the following 
correlation function: 

(WIPW = f(^) (114) 



with t > t'; this becomes, at long times, a scaling function of 

x 



T 2 In t' 

and 



(115) 



a = lnt/lnt' > 1. (116) 

We have computed the scaling function F[X,a] in Appendix C. For simplicity we have set 2g = 1. In Laplace 
transform variables we have 

F]p, a] = J dXe- pX F {X; a) = ^ - ^tanh 2 (117) 
•J 2 — — —coth 2 [^P] [8 + 3p - 16 cosh(^/p) + sinh(^/p) 



a 2 p 3 sinh (^/p) \ 2 

+ (8 + 5p) cosh(2Vp) - 12Vpsinh(2Vp)] 

coth 
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\±-L (1 - VPCoth(Vp)) ( ^=tanh (^) - 1 



a 2 p 2 y/p - VVP V 2 

U> (1 - Vpcoth(Vp)) 2 



a 2 p3sinh 2 

Because the change of correlations between two points is caused only by passing of domain walls through one of the 
end points, the two-time correlation function decays at large x to the square of the autocorrelation function (105), 
i.e \im p _ >0 pF[p, a] — (4a — l) 2 /(9cn 4 ). The spatial decay to this constant value is determined by the closest poles 
to the imaginary axis in the complex p-plane. This yields exponential decay with a characteristic length which is the 
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maximum of £(t') = (T' /tt) 2 and £(t)/4 = (T/2tt) 2 . In the regime a ~ 1, corresponding to lnt w hit', we have the 
following expansion to order 0(a — 1): 



/ dXe- pX F(X;a) = - - (a - 1)4(1 
7o P P 



VP 



sin 



(118) 



Note that in experiments one could in principle measure the Fourier transform which is related to the cross- 
correlations of the scattering speckle patterns at two different times. It reads: 



+oo 



hit . 



]T (S (t)S x (t))(S (t>)S x (t>))e^ = {T\nt') 2 H(q = Q{T\nt') 2 ,a=^) 

X — — OQ 

H(q, a) = 21l(F\p, a] - (4a - l) 2 /(9a 4 p))\ p=lq ) q ? 



(119) 
(120) 



where 1Z denotes the real part and the 1/p part has been substracted to get rid of the 5(q) part. The value of the 
Fourier transform as q — > is: 



H(q = 0,a) 



6 + 40 a - 59 a 2 - 20 a A + 45 a 4 
135 a 4 



(121) 



which is plotted in Fig 5. The Fourier transform is plotted in Fig 6 for several values of a. Note the maximum at 
q > which develops for large a and which is related to the nonmonotonic behaviour of the correlation as a function 
of x. The fact that the above correlation indeed reaches its limit by below can be seen as follows for large a. For 
x <C T 2 , and x belong to the same domain and thus the above two points two time correlation is approximately 
equal to < S[)(t')S x (t') >. This however decays exponentially with x/T' 2 (and thus exponentially in a 2 if one chooses 
ai/F 2 fixed but very small) while the asymptotic value (105) decays only algebraically. 



FIG. 5. Value at q — of the Fourier transform H(q = 0, a) (y axis) as a function of a = hit/ hit' (x axis) 
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FIG. 6. Fourier transform H(q,a) (y axis) of the scaling function for the two point two time spin glass correlation as a 
function of q (x axis) for four different values of a = 1.25, 2., 5., 20. 

VIII. RARE EVENTS, TRUNCATED CORRELATIONS AND RESPONSE TO A FIELD 

In this Section we compute time dependent thermal (truncated) correlations as well as the response to a uniform 
field applied at time t w . For this one needs to go beyond the effective dynamics which just places each wall at a 
specific location at each time. Indeed, in the effective dynamics, the local magnetization (i.e the thermal average) 
(S x (t)) at a given point x is given by the orientation of the renormalized bond containing the point x at scale T, and 
is thus either +1 or —1. Thus truncated correlations are zero to leading order and to estimate them one needs to 
consider the rare events in which a domain wall can be found with substantial probabilities at two different positions. 
Such events occur with probability 1/r and the two positions of the domain wall when they do occur are typically 
separated by distance of order T 2 . For example, for the single point Edwards- Anderson order parameter these lead 
to corrections to the zero temperature value of unity of order 

l-«S,(i)»2oc_L (122) 

In this Section to simplify the notation somewhat, we set 2g = 1. 

A. Description of the important rare events 

The rare events that are important for the RFIM turn out to be the same (as far as the energy landscape is 
concerned) as the ones that we considered in our previous study of the aging properties of the Sinai model [28] , with 
a slightly different physical interpretation and different observables to be computed. There are two types [42] of such 
rare events that occur with probability 1/r, denoted (a) and (c) in [28]. We now describe them. 

Events (a) concern bonds which contain two almost degenerate extrema (see figure 7). In the following, we will need 
the probability D^ r , (x) that two points and x > belong to a renormalized bond at T which has two degenerate 
extrema separated by a barrier smaller than T' such that the two points are both located between the two degenerate 
extrema. This is, using the calculations of [28]: 
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A x A 

FIG. 7. Rare events (a) and (c) discussed in the text 



4t<M = 4/ Pr(l-y)(y-x) [ dTo^f^) (123) 

1 Jl>y,y>x JO 1 1 

where the functions f(Y) and G(F) have been defined in (E5-8) and (174) of [28]. 

In such an event (a), the corresponding domain wall (A or B) will fluctuate thermally between the two extrema (see 
figure 7). At equilibrium, the thermal probabilities of finding the domain wall in the two positions are p and (1 — p) 
with p = 1/(1 + er u l T ) and u representing the energy difference between the two minima. The random variable u is 
distributed uniformly around u — 0. In the following, we will need 

M r ) = W T^ = ir4»^^^ = r^ F IM_ (125) 

the factor 1/2 arises as the integral is only [43] over u > 0. 

Events (c) correspond to bonds about to be decimated, 1 
segment [0,x], x > belongs at scale T to a bond of barrier V (i.e. £ = 0) 

9 f 9 f+°° rp 

D t\x) = f2 J MC = 0,l)(l-x) = f J o dX XP*( V = 0, A + — ) ( 1 2(i) 



Events (c) correspond to bonds about to be decimated, with a barrier r + e w V. The probability D^\x) that the 



T 2 ' 



f^r 1 ^-^ (127) 
fe=i 



For the events (c) associated with a barrier T + e, the probability p c that the two corresponding domain walls have 
not yet annihilated at T is given by p c — exp(— e~ c / T ). In the following, we will need 



d n (T) = (4p c (l-p c ))" = T4 n / —e- n m-e- z ) n =T2 2n Y^C k n _ 1 Ml + -—) (128) 

Jo z k + n 

where we have used that the distribution in e is uniform around e = 0. 

B. two point truncated equal time correlations 

1. N on- equilibrium behavior 
We consider the various moments of the truncated equal time correlation function defined as: 
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C n (x,t) = [(S (t)S x (t)) - (S Q (t)){S x (t))] n \ (129) 

C n (x,t) is the sum of two contributions. The first contribution C^(x,t) originates from events (a) described above, 
where we have (So(t)S x (t)) — +1 while (So(t)) — (S x (t)} = ±(1 — 2p), because the domain wall fluctuates between 
the two extrema (see figure 7). We thus obtain the first contribution as 

Ci a \x,t) = c n (T)D^ r (x) (130) 

with r = Tint. 

The second contribution Cn\x,t) originates from events (c) described above, where we have (So(t)S x (t)) = +1 
while (So(t)) = (S x (t)) = ±(1 — 2p c ). The second contribution is thus given by 

C ( n c \x,t) =d n (T)D c r (x) (131) 

The final result for the moments of the truncated equal time correlations is thus: 

c^, t) .^^m^hm.-^ (132, 

fe=i 

where c n (T) and d n (T) are given in (125,128). Note that to this order in 1/r the n dependence is only contained in 
the prefactor and in particular the correlation length £th(t) — (Tint) 2 /tt 2 extracted from (132) does not depend on 
n and is equal to the result (82). 

For the spin glass, the above formula gives the even moments C2 n {x,t). 

For the following sections, we will need: 

x=+oo 

Ci(t)= C 1 (x,t) = (- + -ln2)T 2 lnt (133) 



2. Equilibrium truncated correlation 



To compute the equilibrium truncated correlations the method is very similar. One must stop the RG at scale 
Tj = 4J and consider again events (a) (which give the same contribution as above replacing T by Tj) and events 
(c) but with a different interpretation and result since there are now at equilibrium. Indeed in the renormalized 
landscape at scale Tj the only thermal fluctuations (apart from the (a) events) come from barriers T = Tj + e w Fj. 
The barriers much larger than Tj are occupied by a pair of domains with probability almost 1 while the barriers well 
below Tj (which have been decimated at previous stages) are occupied with probability almost 0. The barriers with 
r = Tj + e w Tj are occupied with probability p = 1/(1 + e~ e / T ). Thus we now have the equilibrium truncated 
correlations: 

C^(x) = c n (T)(Dl rXj (x) + D c Tj {x)) = -c n (T) £ 1 [ 2 J > e~ W (134) 

fc=i 

Here, as above, the correlation length £f? = 16 J 2 /tt 2 extracted from (134) to this order in Yj does not depend on 
n. Since it was argued in [22] that the correlation lengths of the Cn 9) (x) generically depend on n, our results suggest 
that here this dependence is subleading in Tj. Our result for the correlation length of c[ eq \t) coincides with the 
result of [21] and with the result (83) with 2g — 1 [44] However the detailed form of the functions Cn q \x) obtained 
here explicitly depends on n. 



3. Approach to equilibrium for truncated correlations 

Since the equal time result for r < Tj = 4 J and the equilibrium result differ only by substituting d n (T) by c„(T) 
in the (c) events, there should be a non trivial crossover near r = Tj towards equilibrium controlled by events (c), 
which we now analyze. 
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Let us consider a bond with barrier F — Tj + e. When r = Tint is close to Tj this bond can be either occupied 
(with probability p(t)) or empty (with probability 1 — p(t)). One has: 

^ = I(l-p(t))-Ip(t) (135) 

at T\ T2 

where n = e Tj / T is the inverse rate of creation of a pair of domain walls (which immediately migrate to the endpoints 
of the bond) and T2 = e^ r J+e )/ T is the inverse rate of annihilation of the pair of domain walls located at the endpoints 
of the bond. Thus, substituting t = e r / T one finds: 

* t >-i+^ + i^-*-< i+ ^£> <>*> 

where lnt eq = Tj/T and one has t/t eq = e~( r - , ~ r ) / " r . Integrating over e, one obtains that the crossover to equilibrium 
for t <~ t eg is described by: 

C n (x 1 t)^Ci a \x)+C^(x) (137) 
C^(x)=c n (T)D^ rj (x) (138) 
C^{x) = e n {T,t)D<f(x) (139) 



with: 



f+00 7 n—\ 

e n (T, t) = 4p(t)(l - p(i)) = 4"T / dz- — (1 - e -(i+*)t/*.«)"(i + ^-(l+zjt/t.,^ (140) 

Jo (1 + z ) ™ 

n n 

= 4 " T E E(- 1 ) fcC " C ^ r [ n + rf^( n + P' 1 + P - n ' ( fc + P) — )e" (fc+P) ^ (141) 
fe=o P =o te « 

where lni eg = AJ/T. This expression crosses over from e n (T,t/t eq <!)-* d n (T) and e n (T,t/t eq » 1) — > c„(T). 



C. Two-point two-time truncated correlations 

We now consider the truncated two-point two-time correlations 

C n (x,t,t w ) = ((S (t)S x (t w )} - (S*oW)(^(^)»" (142) 

The calculation is very similar to the equal time truncated correlations. We first consider the events of type (a), 
where we have now to keep track of the barrier T between the two almost degenerate extrema. There is a non- 
vanishing contribution if the barrier To is smaller than r„, but bigger than T = Tln(t — t w ), so that equilibration 
cannot take place between t w and t. In that case we have (So{t w )S x (t)) = +1 while {So{t w )) = (S x (t)} = ±(1 — 2p), 
where p = 1/(1 + e~"/ T ) as introduced above in the description of events of type (a). These events lead to the 
contribution 

C^(x,t,t w )=c n (T)± j P r (l - y)(y - x) f~ dT ±f(^) (143) 

Jl>y,y>x Jt i o 

= c n (T)(D a TXw {x)-D a Tt {x)) (144) 

where c n (T) and -Dp r (x) are given above. Note that the contribution of these events vanishes when T becomes 
equal to r„,. 

We now consider events of type (c), which give different contributions and must be examined separately in the 
scaling regime T = Tln(t — t w ) < T\nt w = T w and the scaling regime t — t w <~ t w , i.e T = T W . 

We first consider T < T w . In that regime the events of type (c) should be considered at scale T w where we 
have (So(t w )S x (t)} = +1 while (So(t w )) — (S x (t)) = ±(1 — 2p) which, together with the (a) events, gives the total 
contribution 

C n (x,t,t w ) = Ci a Hx,t,t w )+d n (T)D c r Jx) (145) 



2G 



this formula holds for f = Tln(t - t w ) <Tlnt w =T W . 

In the regime t — t w <~ t w i.e T = T w , the events (c) also start equilibrating, which we now study. Let p c (t w ) = 
exp(e _e / T ) the probability that the domain walls separated by the barrier r„, + e have not yet annihilated at t w . Let 
p c (t) — p c (t w ) exp(— t-^e~ f -/ T ) = cxp(— j-e~ f -/ T ) the probability that in addition they have also not yet annihilated 
at t. One has (with x and belonging to the bond being decimated): 

< S x (t w ) >= 1 - 2p c (t w ) (146) 
<S x (t)>=l-2p c (t) (147) 

< S {t)S x (t w ) >= 1 - 2p c (t w ) + 2p c (t) (148) 

< S (t)S x (t w ) >-< S (t) >< S x (t w ) >= 4 Pc (t)(l -p c (t w )) (149) 

Thus one gets, in the regime t — t w ~ t w the total contribution 

C n (x,t,t w ) =d n (T,t,t w )D°Jx) (150) 

d n (T, t, t w )) = (4p c (t)(l - p c {t w )Y = T 4" / -e-"^ z (l - e~ z ) n (151) 

Jo z 

In the limit t >• t w these truncated correlations decay, for fixed x/(T\nt w ) 2 , as C n (x, t, t w ) ~ (t w /t) n . They decay 
to zero as there are no other contributions for later time t. 
In the following we will need 

x— +oo 

Ci(M«0= I] (152) 

a:— -co 

= 2! T>( „_(M£zW)! )+ » T!ln21I , t „ (153) 

45 (mttu)^ 45 mt w 

= | T2lnl . ln(1 + k) t - 1 ^- 1 P") 

= ^mt. for !£(i_M>i ,155) 

45 i lni„, 



D. Response to an applied field 

In order to compare with typical aging experiments, we will consider the following two histories for the system and 
compare them: 

(i) Apply H > starting from t = : the magnetization per spin m(t) will then grow in time as computed in (60) 
up to time t eq ~ e 4J ' T , where m eq is reached. 

(ii) Keep H = between t — and t Wl and then apply H > for t > t w : in this case the magnetization per spin 
m{t,t w ) remains up to time t Wl and then grows to again reach m eq in the large time limit. 

We now estimate m(t,t w ) in the case (ii) in the "small applied field regime", where H ~ 1/T^,. It is convenient to 
define 



and introduce the ratio between 



and r,„ = T In t,. 



r^, In t w 

We discuss separately the three regimes 0<&<1,&~1 and a > 1. 



i=t-t w (156) 



r = Tlnt = Tln(i-i w ) (157) 



f _In(*-*»). (158) 



27 



1. Response at early times < a < 1 from degenerate wells 



We first study the scaling limit of small H and large r„, with HT^ fixed and < a = < 1 fixed. In this regime 
the dominant contributions come from bonds with near degenerate extrema, i.e the rare events of type (a) described 
in (VIII A). 

Let us consider an ascending bond with a secondary minimum separated by a distance y and a barrier T at a 
potential u > above the minimum. Using the results and notations of [28] we find that the probability that an 
ascending bond of this type at T w (before the field is turned on) is: 

r r „ (y, T )dT dy = 29(T W - r )P r „ (0, .) * y P r „ (0, .)dT dy (159) 

Just before T w , there is thermal equilibrium between the two extrema and thus the probability that the primary 
cxtrcmum is occupied by an A domain wall is p eq {u) = 1/(1 + e~ u ' T ). The thermally averaged total magnetization 
of the segment y is then < M >= y(l — 2p eq {u)). One now turns on the field at T w and the thermally averaged 
magnetization of the segment remains the same (up to negligible probability) until the time r = Tq when a new 
equilibrium is attained (here we can neglect the equilibration time scale). Just after T — Tq the occupation probability 
of the left minimum is now p eq (u') = 1/(1 + e~" / T ) where u' = u — 2Hy. The new magnetization is < M >'= 
y(l — 2p eq (u — 2Hy)). Similarly, the contribution of a descending bond (with y, Tq) is given by the same formula with 
u < 0. Finally the contribution of degeneracy of hills — fluctuations of the B domains — yields an overall factor of 
two. The total contribution of all these events to the magnetization per spin is: 

^ rT r+oo r+oo 

m{t,t w )=m^ a \t,t w ) = 2- T dV dy du r Fw (y, V )y(l - 2p eq {u - 2Hy) + 1 - 2p eq (-u - 2Hy)) (160) 

1 w J0 JO JO 

This gives: 

^S(HT 2 ) — — f dYY 2 G(Y) (162) 



r r 

-45^rI-45 (HrJ f^ (163) 

In this regime, the magnetization as a scaling function of HT^ is thus exactly linear. It can be shown that non 
linear response in HT^ couples to the curvature of the distribution of difference of potential near zero and is of higher 
order in 1/T W . 



2. Response at times a ~ 1 

When t — t w is of order t w a second effect adds to the one computed above. It corresponds to the events of type 
(c) described in Section VIII A where the barrier of a bond at T w is equal to T w + e where e = O(l) (of arbitrary 
sign). In the absence of the field the pair of domains at the endpoints of the bond have not yet annihilated at time 
t w , with a probability p c (t w ) = exp(— e _£ / T ). When adding the field at t w the barriers suddenly either increases (for 
descending bonds) or decreases (ascending bonds) by 2HI where I is the length of the bond. For t > t w (and such 
that t — t w ~ 0(t w )) the probability p c (t) that the domain has not yet annihilated depends on H and is: 

Pc(t) - Pc (t w ) exp(-(^) e (-^>/ T ) (164) 

tw 

for descending and ascending bonds respectively. The events (c) thus result in a difference in magnetization compared 
to the zero field case equal to: 

m (c) (Miu)= 2 [°° dl / + "d e p*(0,0?exp(-e- £ / T )(cxp(-^^ e T(- e - 2 ^)_ cxp( _Ll^ e ^(-+2^))) (165) 

I w JO J -co t w t w 

2T r°° t +(t-t ) e 2HXV ™/ T 

-t.L w '"- t ' i " i| ,:Ut-"^ 1 (w6) 
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The total magnetization is now: 

32 

m(t,t w ) = m {a) {t,t w )+m {c) (t,t w ) = —HT\nt w +m {c) {t,t w ) (167) 

Note that in the present regime the magnetization as a scaling function of HT^ is complicated and non linear. 
In the limit where HT^ is small, one has 

m^(t,t w ) = ±m*J=± r d XP* (ll = 0, A) A 2 + 0{{HTlf) = ^HrJ—^ + 0((HTlf) (168) 
1 to t J Q 45 t 

where we have used J™ d\P*(r] = 0, A) A 2 = 7/180. 

Although the above function (166) is complicated, at the special time such that t — t w = t w it takes the simple 
value: 



4 f°° 7 

m^(2t w ,t w ) = —HTl / d\P*( v = 0, A)A 2 = —HT V 

I to Jo 45 



(169) 



3. Response at times a — a > 1 

For time differences t — t w ^>^> t w , corresponding to a — a > 1, the response of the RFIM chain to an applied 
field will be dominated by the effective dynamics described by the RSRG procedure. When the field H > is turned 
on, the descending bonds with (F,l) become (F + 2Hl,l) and the ascending bonds become (F — 2Hl,l) except if 
< F < T w + 2HI, since in this case they must be immediately decimated. Technically, from the point of view 
of the landscape, it is more convenient to symmetrize the initial condition at T = T w which amounts to artificially 
reintroduce the descending bonds r„, — 2H I < F <T W (these bonds being redecimated immediately do not introduce 
any errors for a > 1). This correspond to the following initial distributions at T = T w : 

/•OO 

P+ {F,l) = Pt w (F, I) - 2HW F Pr w (F,l) + 2H{Pv w * FA Pr w n (IPrjOJ)) - AHPr w (F,l) / dl'l'Pr w (0,l')) (170) 

Jo 

POO 

Pf (F,l) = P Tw (F,l) + 2HW F P rw (F,l) - 2H(P Tw * F> , P Tw *, (lP Tw (0,l)) + 4HP Tw (F,l) / dl'l 1 P T J0,1')) (171) 

Jo 

We now check that the magnetisation corresponding to this initial condition is the one at the end of the a — 1 regime 
(i.e the t — ► +oo limit of 168). It is, to first order in H 



m e// (f = T w , T w ) = J ° ^ tV^TTTT, , ( 172 ) 



j +oo dll^dC [P?J(,i)-PfJt,i) 
f +co dllf °°dC[Pl(t,i) + P r -JC,i)] 

Q TT POO 1 A 

dllP Tw {T w ,l) = -HT w (173) 

1 w JO ^ 

We write to linear order Pp(F, 1) = Pf,(F, I) + 2HQt;(F, I) with the initial condition in Laplace transform variables 

tit r =: r.^ I 

Q+JCP) = -d P (Ur a (p)u Tw (p)e-^W) - U^{0)UrApfC^ Ur - {p) + 2U^jO)U r Jp)e-^^ (174) 
Q^JCp) = -Q+JC,p) (175) 

To compute the magnetization we are interested only in P^(F, l) — P7(F, I) — 2HQ^{F,l) where Qf.{F,l) — Qf(F,l) — 
Qp (F, I) satisifes the linearized RG equation: 

(0 t - d ( )Qr(C,p) = -Q t (0,p)P t (.,p) * c P t (.,p) + 2P f (0,p)Q f (.,p) * c P f (.,p) + 2Q f (0, 0)P f (C,p) (176) 
The solution has therefore the form: 

Q t (C,p) = (Af(p) + CB t (p))e-^r(p) (177) 
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where the coefficients Af,(p) and -Bf (p) satisfy the RG equations 

d t A t (p) = -u t (p)A t (p) + B t (p) + 2U t {p)A f (0) 
d t B f (p) = -u t {p)B f {p) 

with initial condition at Y = Y w 

A Tw (p) = -2d p (U Tw (p)u rw (p)) + W Tw (p)U^ (0) 
B r Jp) = 2U r Jp)u r Jp)u' r Jp) - 2U r J P ) 2 U^J P ) 



The solutions are 



B t (p) = Br a (p)- 



sinh(T wy /p) 



sinh(r^) 



A f (p) = (A r Jp) + (T - r w )B r Jp))- 



sinh(r„,^/p) 



In the limit f 3> r„,, we have 



B f (p) ~ - 



sinh(f y/p) 

Vp 



2(t-T w )- 



sinh(rvp) 



A t {p) * -T 



sinh(rvp) 



sinh(r^p) 

and recover the first order linearisation in 5 of the biased fixed point solutions (37) 

P r ± (C,p) = U±(p)e-rtM = Pt(C,p)(1 T 5T)(1 ± «C) = ft(C,p) + 2(«C - ^ r ) sinh ^ } ' 
The magnetisation at first order in H is given by 



-Cv^coth(rvp) 



m e//(f) 



r 2 7o 







> O 7^ 

dCQf(C.0 = -7 



A t {p) | B f (p) 
«f(p) «f(^) 



Ip=o 



2 _ 16 /rv 

3 45 V f 



It grows from 



m e//(r = i\ 



14 

45 J 



to the asymptotic regime for large f >f, 



m eff (t > T w ) = -HT 



that corresponds to the behavior of the magnetisation of the biased case at first order in 5 (60). 
We now summarize our results for the magnetisation in the regimes a < 1, a ~ 1 and a > 1, 
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m(r,r w ) = — m^a 3 



for 

t tq. 



a= — <l 



32 14 
m(T,T w ) = -HT w + -HT w ; 



i(r, i\„) = ht 



2 16 

3 ~ 45a 3 



45^ 



for a = — — ~ 1 

A in 



for a = — > 1 
J- to 



(178) 
(179) 

(180) 
(181) 



(182) 
(183) 

(184) 
(185) 

(186) 

(187) 

(188) 
(189) 

(190) 
(191) 

(192) 
(193) 
(194) 
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E. Fluctuation-dissipation violation ratio 



Having computed truncated correlations and the response to an applied field, wc now discuss the fluctuation- 
dissipation violation ratio, a measure of the nonequilibrium behaviour of the system. For two observables A and B, 
the fluctuation-dissipation-theorem (FDT) violation ratio X is defined as [6,45] 

T R A ,B(t,t w ) = X A>B {t,t w ) d tw C A , B {t,t w ) (195) 

where CA,B(t,t w ) represents the truncated correlation 

CaAWw) =< Mt)B{t w ) >-< A(t) >< B{t w ) > (196) 

and B,A,B(t, t w ) represents the response in the observable A at time t to a field Hb linearly coupled to the observable 
B in the hamiltonian through a term of the form —HbB 

RA,B{t,t w ) = ——jj—\ Hb=0 (197) 
otl B {t w ) 

Here we have computed the magnetization resulting from an uniform magnetic field H, so that the observables A 
and B are given by A = (1/L) an d B — J2i=i &i respectively. From the magnetization 

t 



m(t,t w ) = - 2^ < Si(t) >= < S (t) > = H I duR{t,u) (198) 



and the truncated correlation 

L L 

L 



^ L L 

d(t,t w ) = jJ>2J2( < ^W^(^) > - < Si(t) >< Sj{t w ) >) = J2 < S (t)S x (t w ) >-< S (t) >< S x (t w ) > (199) 

i— 1 j — 1 x 

one obtains the fluctuation-dissipation ratio X(t 7 t w ) as 

Note that we have used the infinite size limit to replace translational averages by disorder averages. 

We have computed C\(t, t w ) (153-155) and m(t,t w ) (192-194) in the three regimes < & < 1, a ~ 1 and a > 1, 
and thus we get the following expression for the fluctuation dissipation violation ratio at large times t,t w ^> 1: 

X(t,t w ) = l for < a = ln(t ~ tw) < 1 (201) 

In t w 

X{t 1 t w ) = t -^ for — — a fixed number (a = 1) (202) 

t t w 

x(t, '" ) = t^: (1+ yi?r ) fora= ^- > 1 (203) 

The behaviour of X upon increasing the time difference t — t w is as follows. First of course one expects, when 
t — t w is a fixed number, a non universal equilibrium regime (not studied here) where time translational invariance 
and FDT is obeyed. Increasing the time difference as t — t w ~ in the regime < a < 1, we find that the FDT 
theorem is still obeyed to leading order, which is quite interesting since time translational invariance does not hold in 
that regime. Next, in the regime a = 1, t/t w fixed, X becomes a non trivial scaling function of (t — t w )/t w , which 
interpolates between X = 2 for (t — t w )/t w — > and X = 1 for (t — t w )/t w — ► +oo. In order to match the value X = 1 
there is thus a non trivial crossover regime between the end of the quasi equilibrium regime a < 1 and the beginning 
of the a = 1 regime. This crossover occurs for (t — t w ) <~ lnt^, and is given by 

X(t,t w ) = F 1 {{t-t w )lnt w /t w ) (204) 

F l(y) = ^±^ (205) 
v ' 7y + 96 v ; 
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Finally, for very separated times, in the regime a > 1, we find that X grows towards +00. This occurs again after 
a crossover between the end of the a = 1 and the a > 1 regime which occurs on time scale t ~ t w In t w , where X is 
given by 

X(t,t w ) = l + ^-J— (206) 
' m i w 

which matches both the required limits. 

We can now compare with the mean field models [49]. As in mean field [6,45], we find here an aging regime where 
X is non trivial, and for t/t w fixed it is a function of this scaling variable. On the other hand, contrarily to mean field 
models, the ratio X here is never a function of only Ci(t, t w ) 1 in the regime d = 1 because of the extra power of hit in 
Ci(t, t w ), and in general because of the presence of both scales (t, t w ) and (lni, \nt w ). In addition, X is found here to 
be nonmonotonous and the values taken by X are not within the interval [0, 1]. In particular X tends to +00 in the 
asymptotic regime t ~ i™ with a > 1 since truncated correlations become very small compared to the response in that 
regime. Since the ratio X has an interpretation in some contexts as an inverse effective temperature X = 1/T e ff [45], 
one would find here that T e ff — ► at large time separations, in contrast to the result that T e ff — > +00 in mean field 
models. Although this might appear surprising at first sight, one should remember that in finite dimensions many of 
the properties of the RFIM are controlled — in the the renormalization group sense — by a zero temperature fixed 
point; this includes the intermediate regime of length scales studied here for a system — the ID RFIM — with no 
phase transition. 

It is also interesting to compare our results with the fluctuation-dissipation ratio X in pure systems presenting 
domain growth, typically ferromagnets below their critical temperature T < T c [7-9] or at T = T c . The ratio X is 
usually computed by considering the local observables A = B = So but to compare with the present study we will 
translate these results for a spatially uniform applied field. 

Let us briefly recall why X = in the large time scaling regime for T <T C (X decays to for large t w ). When initial 
conditions (high temperature T > T c ) before the quench are chosen to have zero magnetization (J2 X < S x (t) >= 0) 
the spin autocorrelation coincides with the truncated correlation J2 X < S x (t) >= and takes a simple scaling form 
[1] 

Cr e (t,t w ) =^2< S (t)S x (t w ) >^ L{t w ) d h{^-) (207) 

where L(t) is the typical size of domains at time tanrf the dimension of space. On the other hand, the magnetization 
when a uniform field has been applied between t w and t behaves as [8-10] 

M auto (t, t w ) ~ L(t w ) d - a f 2 (^) (208) 

with a = 1 for Ising order parameter [7-9] and a = d — 2 (d > 2) for O(N) model [9]. Note the reduction of M 
with respect to C\ by a factor l/L(t w ) in ferromagnets, which immediately implies X = in the scaling regime 
L(t w ) ~ L(t) (X = as soon as a > 0). As is usually argued, the origin of this factor can be seen by considering 
the system at t = t w and focusing on the immediate response to a small applied pulse field [10] (here in d = 1 for 
simplicity): each interface responds by a factor 0(V) and since they occupy only a portion 1/ L(t w ) of the system this 
yield a total response 1/L(t w ). By contrast, note that exactly at criticality, T = T c , there is a non trivial X = X(t/t w ) 
which appears as a dimensionless amplitude ratio, by a different mechanism [9] . 

For the RFIM case, the full correlation function has the same scaling form as (207) (as can be seen by generalizing 
the result for the autocorrelation computed in (105)) but the truncated does not since we have obtained (153-155) 

Ci{t,t w ) = lnt w M H *~ tw) ) for < ln( *~ U < 1 (209) 
in t w in t w 

= lnW 2 (y) for ^ ~ ^ > l (210) 

whereas for the RFIM magnetization we have found (192-194) 

m(t,T w ) = lnt w M Ht ~, tw) ) for & = ^- < 1 (211) 

in i w 1 w 
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m{t,r w ) = ]nt w ip 2 ( t -f) for a = - 1 ( 212 ) 

m{t,T w ) = \nt w ^{^-) for d = > 1 (213) 

These expressions are rather different from the ferromagnet (207), first because both t and Int appear in these scaling 
forms. The only domain growth-like non trivial scaling regime with L(t) = T\n 2 t occurs for a > 1. In this regime 
the magnetisation has a form similar to (208) with a reduction factor 1/ -J L(t w ) (rather than 1/L(t w ) in the pure 
system). Using the discussion of Section VIII D 3 one understands that the origin is quite different from the pure 
case. In the immediate response to an applied field (173) only a small fraction of domains 1/T W ~ 1/ \jL(t w ) respond 
but their response is very large since the full domain, of length ~ L(t w ) flips. The truncated correlation on the other 
hand is very small and does not take the scaling form (207). This yield a value X = +oo in this regime. 

Similarly, the origin of the non trivial value of X in the regime t~t w (a=l)is very different from the pure 
case. Both the response and correlations originiatc from rare events and take the form ln * L(t w )f(t/t w ) where now 
the factor -^j— is the probability of the rare event and f(t/t w ) its contribution to activated dynamics. Since they are 
of the same order this yield a non trivial X. 



IX. PERSISTENCE PROPERTIES OF THE RFIM 



We now turn to a study of the persistence properties of the random field Ising model which have received substantial 
attention for other systems evolving towards equilibrium. Two of the primary properties of interest in this context are 
the time decay of the probability that a spin has never flipped up to time t and the time decay of the probability that 
a domain has survived up to time t. The results for the single spin persistence for the RFIM (sections IX A and IX B 
below) are also valid for the spin glass. The large time limit of these quantities can be computed from the asymptotic 
full state (see fig. 3) on which we focus below. 



A. Persistence of a single spin 

In zero applied field the probability n(t) that a given spin at x — has never nipped up to time t in a single run, 
is equal to the probability that neither the nearest domain wall on one side, nor the nearest (opposite type) domain 
wall on the other side have crossed x = 0. 

In [28], we have found that the probability Ili(t) that a given Sinai particle does not cross its starting point up to 
time t decays as III (t) ~ l{t)~ Bl with 6\ = 1/2. We thus obtain that II (i) in the zero field RFIM decays as 

U(t) ~ I(t)- with = 20i = 1 (214) 

This should be compared with the result [46] for the pure Ising model, where = | (corresponding there to the charac- 
teristic length l{t) ~ \/t). This out of equilibrium behavior holds up to time t = t eq corresponding to renormalisation 
scale r = V j at which equilibrium is reached. 

In the presence of an applied field H > 0, we can use our previous result [28] for the biased Sinai diffusion. We 
have found that the probabilities 11+ (t) (resp. 117/ (t)) that a given Sinai particle remains on the right (resp. left) of 
its starting point up to time t, in the case of a drift in the (+) direction, behave as 

nfw-r^r (215) 
n r(')*^r7 ( 216 ) 

This gives the the probabilities II * (t) that a given spin in the RFIM keeps the value (±) up to time t, which in the 
limit of large ST In t » 1 behave as 

n+(t) « ^ (217) 

U-(t) w ^t- 2HT 'a (218) 
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B. Persistence of the local time-averaged magnetization 



In addition to the persistence of a single spin, one can also obtain the statistics of the flips of the thermal average 
of the local magnetization, i.e (S x (t)}, at a given site x. As explained in [28] in the context of the Sinai model, 
quantities averaged over many runs behave very differently than quantities for a single run; in particular while the 
spin of interest in one run may flip many times, if the same spin is examined in many runs with statistically similar 
initial conditions, the average over the runs at a given time into the runs will flip far less often. 

The local magnetization will successively be equal to I and —I with only very small probability, at large times, of 
being be equal to an intermadiate value. The sequence of flips is given by the sequence of changes of orientation of 
a bond during the renormalisation procedure extensively studied in [28]. From that analysis, we obtain the following 
results for the RFIM. 

In zero field H = 0, we denote by k the number of changes of sign of (S x (t)) at a given point x between and t. 
The distribution of the rescaled variable 



is characterized by the asymptotic decay 



lnT ln(Tlnt) 



(219) 



where the generalized exponent 9(k) is 



(?(«) = -In 



Prob( K ) ~ ~i(ty 











2k 
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K 2 + -. 



(220) 



(221) 



The exponent 0(k) is a positive convex function : it decays from 6(k = 0) = 3 ^ to 0{\) = 0, and then grows again 
for k > 1/3. This implies in particular that 

k 1 



K = 



ln(T hit) 



with probability 1 at large time 



(222) 



All of the moments of k will be dominated by the typical behavior; i.e. {n m ) = 3 m for all m. The full dependence 
on k of the 8(k) function describes the tails of the probability distribution of k, i.e., the large deviations. Note also 
that the probability that (S x (t)) doesn't flip up to time t decays as 1{t)~ e with exponent 



6 = 6(k = 0) 



V5 



(223) 



which is significantly smaller than the corresponding decay exponent 8 = 1 in (214) for a single spin. 

Since the renormalisation procedure has to be stopped at T = Tj at which the equilibrium is reached, and that at 
later times no more changes occur in the local magnetization we obtain that the total number of flips is: 



hot -» ^ ln(4 J) 



(224) 



the decay of the tails of the probability distribution of k = ktot/ In (4 J) being described by the same function 0(n) as 
above in terms of the length Lim- 

Another result from [28] is the characterisation of the full sequence of the times T\ = Tlnti, ...Tk = Tlnt^ where 
the local magnetization (S x (t)) flips. The sequence of scales {1^} is a multiplicative Markovian process constructed 
with the simple rule Tk+i = ctkTk where {ctk} are independent identically distributed random variables of probability 
distribution p(a) 



p(a) = 



1 



a X + — A_ 



(a 



— a 



with 



A± = 



3± V5 



(225) 



As a consequence, Tk — ak-iotk-2 • • • is simply the product of random variables, so that we obtain using the 
central limit theorem that 
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fe^oo \ k J 



(226) 



and we thus recover that the number k of changes of (S x (t)) grows as lnT = In hit and that the rescaled variable 
k = ln ( T fc lnt ) is equal to 1/3 with probability 1, as in (222). 

These results can be extended to the RFIM in the presence of an applied field H > 0. The total number of flips in 
this case saturates to a finite value given by a scaling function of H and J identical to the one given in Section IV E 
in [28] for the total number of returns to the origin in the Sinai model. 

C. Domain persistence 

Persistence can also be defined for larger scale patterns. In the RFIM the simplest pattern (beyond a single spin) 
is a domain. When a domain of e.g. consecutive + spins disappears, the two nearest domains of — spins merge. Thus 
domains can either grow by merging with other domains, or die, and this leads to interesting survival properties, that 
were studied by Krapivsky and Ben Nairn [47] for the pure Ising model. Here we slightly change their definitions of 
the exponents to adapt to the logarithmic characteristic length scale I(t) = 7r ~ (Tint) 2 of the RFIM. 

t=0 




FIG. 8. Temporal evolution of domains in the RFIM. For each surviving domain at time t, m denotes the number of ancestor 
domains in the initial condition at t — 

Of all the domains which still exist at time t, we define R m (t) as the number of domains which were formed out of 
m of the initial domains. This is illustrated in figure 8. Then ^ m R m (t) — N(t), the total number of domains at time 
t, and the fraction of initial domains which have a descendent still "alive" at t is S(t) = J2 m mR m (t) =< m > N(t). 
Following [47], one defines the exponents: 

Nit)^!^)- 1 (227) 
S(t) ~ I(t)-* (228) 
Ri(t) ~~l(ty 5 (229) 

in terms of the mean domain length l(t). Asymptotically, R m (t) has the scaling form: 

iUt)«7(t)*- 2 fl(=^) (230) 

The scaling function behaves, in the pure case, as R(z) ~ z a for small z and as an exponential at large z. One finds 
the exponent relation 5 = v + [y — ip){l + a). 

Let us now study the RFIM using decimation. In the asymptotic state (full state), domains coincide with bonds, 
and when a bond is decimated the domain on this bond dies, and the two neighbors merge into a new domain (the 
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new bond). We thus associate to each bond an auxiliary variable m counting the number of initial domains which are 
ancestors of the domain supported by the bond. Upon decimation of bond (2), the rule for the variables m is simply 

m! = mi + m 3 (231) 

It turns out that this rule coincides with the magnetization rule of the RFTIC [32] and thus one has the scaling 
to - T* with cj) = (1 + \/5)/2 leading to 

i/j = 3 ~ 4 ^ = 0.190983.. (232) 

This should be compared with the value for the pure Ising model, which has only been determined numerically [47]. 
as -0 w 0.252. 

Here it is interesting to note that we have found that tp — 9, i.e. the general exact bound ip < 9 is saturated [29]. 
This bound for tp comes from the fact that a point that has never been crossed by any domain wall up to time t for 
the effective dynamics, necessarily belongs to a domain that has a descendant still living up to time t. In the effective 
dynamics in the "full" renormalized landscape, we have also the reciprocal property : a domain wall surviving between 
t' and t necessarily contains points that have never been crossed by any domain wall between t' and t. Note that the 
equality ip — 9 also appears in the coarsening of domains for the ID T = Landau- Ginzburg 4> A soft-spin Ising model 
[36] for the same reason. The strict inequality ijj < 9 requires the possibility that a domain wall can survive between 
t' and t even if all points inside it at t' get crossed by a domain wall between t' and t (see [29] for examples). 

We now study the probability R\(t) that a domain has survived up to time t without merging with any other 
domain. This requires that the two domain walls A and B living at the boundaries of this domain do not meet 
any other domain wall up to time t. In the asymptotic full state, these conditions imply that three consecutive 
bonds cannot be decimated. Thus the domain in the middle cannot grow and the probability that it survives upon 
decimation thus decays exponentially in T. 

We thus obtain that R\{t) decays exponentially in V (and thus with a nonuniversal power of time determined by 
the initial condition which determines the convergence towards the asymptotic state). It thus corresponds to: 

S = +oo (233) 

As a consequence, the scaling function R(z), not computed here, has an essential singularity at the origin. 



X. FINITE SIZE PROPERTIES 



The real-space renormalization procedure can also be used to study analytically large finite size systems [33,28]. 
We will use extensively the analysis of the Sinai walker with either absorbing or reflecting boundary counditions [28]. 
For the random field Ising model, one may consider several boundary conditions : 

(i) fixed spins at the ends: this corresponds for the Glauber dynamics of the RFIM to domain walls which behave 
as Sinai's walkers with reflecting boundary conditions, 

(ii) free boundary conditions: this corresponds to absorbing boundary conditions for the domain walls. 
We give here some results for (i), whose derivation is detailed in Appendix D and discuss (ii) at the end. 



1. Fixed spins at both ends 

Assume for definitness that spins at both extremities are fixed to the values 5*0 = +1 and Sl = — 1, where L is 
the system size. There is thus an even number of domains in the system, and one can describe its statistics at large 
time using the following generating function. Let us define Ih{k\ T) as probability that the system of size L at time 
t (i.e at scale F = Tint) contains exactly n = 2k + 2 domains, with k = 0, 1, 2... One obtains (see Appendix D) the 
generating function in Laplace transform with respect to the length L as: 



where p = £j. 



' d Le-^(±^ L{ k;r))= 2 gm^Tg) 



(234) 
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In the case of zero applied field H = (S = 0) one easily performs the Laplace inversion and obtains: 

where a = ArcCos^fz E (0, ^) for z G (0, 1). In particular we obtain the average and mean squared total number of 
domains in the system at time t: 

<n>=2<k > L . r +2= Ag-^ + \ + o(L°) (236) 

I z 3 

<n 2 > LX -{<n> LX ) 2 =^ + O(L ) (237) 

with T = Tint. Again, the same quantities at equilibrium arc simply obtained by setting F = 4 J in the above 
formulae. 

In the presence of an applied field H > (5 > 0) one obtains similarly: 

< n >= 2 < k > L r +2 = 4 9 i^^ + 2 1 ~ 7C f h7 + 2 + o(L°) (238) 

sinh 7 sinh 7 

< , 2 > L , r -« n > L . T f = 8gL-^- ( 1 + ^_J^jl) + 0(L o } (239) 

sinh 7 \ smh 7 / 

with 7 = ST Int. 

We have also obtained the probability distribution Fl(M;T) of the total magnetization of the system: 

L i=2fe+2 

M L =J2<Sj>= ( 24 °) 

j=\ i=i 

In Laplace transform with respect to L and M it is: 

[ +OC dLe-^ [ +L dMe-^F L (M;T)^l r ^ + ^ ( _ g " . (241) 
•/ — L 1 - (<? + r)P r (q - r) 

where: 

^ r ^ ~ sinh((5r)(v9+^coth(r v / ^T52) t <s) ^ 242 ^ 

(9) = VP + ° J (243) 

sinh(v5T^L)(VpT^coth(rv^T^) t 5) 

where V = ^g- Note that g + r and g — r are simply the Laplace variables associated to the total rcscaled positive 
length and negative length. 



2. Free boundary conditions 

Free boundary conditions in the RFIM correspond to absorbing boundary conditions for the diffusing domain walls. 
However the study is slightly different from the one carried in [28] because now there are particles (A or B) both at 
bottom and tops. 

The structure of the renormalized landscape at T and the full state near the boundary is now the following (see 
figure 9). There is an absorbing zone of length Iq. Then there is a first bond (barrier Fi length li) of arbitrary sign 
(contrary to Sinai's case where the first bond was always descending). The first particle at the common endpoint of 
the first bond and the second bond and is either a A (descending first bond) or a B (ascending first bond). The RG 
procedure is unchanged in the bulk. The new case is when bond 1 is decimate. Then the absorbing zone becomes of 
length l + h, while the domain wall (A or B) gets out of the system (as the cluster formed by the absorbing zone and 
the previously first bond flips). The bond number 2 keeps the same barrier and length (and becomes the new bond 
number 1). 
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A 

FIG. 9. RG rule near a free boundary as explained in the text. The same rule holds when the first renormalized bond is 
descending by exchanging A and B. 

Let R ± (l) be the probability distribution of the length I of the absorbing zone. It satisfies the RG equation 

d r B%(l) = P£{Q,.)* l I%(.)-B$(l) J dl'P*{0,l') (244) 

In the symmetric case the solution in Laplace transform is: 

i?r(p) = f ^tanh(^) (245) 

Interestingly the shape of the size distribution is the same as for the Sinai particle with absorbing boundaries but 
with a global length rescaling by a factor 1/4. A similat study can be made in the biased case. 

Let us close by noting that the approach to equilibrium will be modified near a free boundary, as compared with the 
bulk (discussed in Section ). Indeed near the free boundary it is possible to create a single domain wall with energy 
cost 2J. Thus, for times such that 2 J < T = Tint < 4 J one must consider different rules for the first renormalized 
bond. 

XI. DISCUSSION AND FUTURE PROSPECTS 

In this paper we have shown how a powerful real space dynamical renormalization group method can be used to 
study the properties of the one dimensional random field Ising model and ID spin glasses in a field. We have recovered 
many known results for the long time equilibrium behavior and obtained some new ones. But the main advantages 
of this method is that it enables the computation and physical understanding of many non-equilibrium features of 
the coarsening process. Although the RSRG method is approximate, it retains all the information needed to obtain 
exact results for universal long time, large distance quantities. As exemplified by the calculations of the response 
and truncated correlations, it has the great advantage over many more formal methods by virtue of providing a 
clear physical interpretation of the origins of the processes (e.g. particular types of rare events) that dominate many 
properties of interest. 

In this last section we address two issues: the potential observability of some of the one-dimensional physics that 
we have studied in detail; and the crucial issue of which features of these one-dimensional systems might apply in 
higher dimensional random systems. 



A. One-dimensional systems 

Although random field Ising ferromagnets are not directly realizable experimentally, other systems in the same 
universality class should be. For random antiferromagnets, an applied magnetic field couples to the antiferromagnetic 



38 



order parameter as a random field, most simply if the local magnetic moments vary randomly; such systems have 
the obvious advantage of the strength of the dominant randomness being readily adjustable. Similarly spin glasses, 
although unfrustrated in zero-field, do exist in quasi-one-dimensional systems. In both of these types of systems, rapid 
quenches could be performed by decreasing the field at low temperatures down to a value of order T <C J at which 
coarsening can take place. In practice, other types of randomness in, for example, antifcrromagnctic or frustrated 
two-leg spin-ladder systems in an applied field may have advantages. 

There is one physical phenomenon of which one must be wary: in many random systems — especially spin-glasses 
[5] — the equilibrium states toward which the system strives depend hypersensitively on macroscopic parameters such 
as the temperature and the magnetic field; this is often referred to — somewhat confusingly — as "chaos" [48]. If this 
were the case here, then changing the temperature or the magnetic field would not correspond simply to speeding up 
the dynamics, but would instead drive the system towards a new equilibrium which might differ on the scales being 
probed from the original one. Indeed, this effect is the origin of much of the most interesting aging phenomena seen in 
three-dimensional spin glasses [5,3]. Fortunately — although less interestingly — this effect does not occur in random 
field chains: although the effective local random fields could change in a non-uniform way with temperature or with 
applied field, this will lead only to smooth modifications of the large scale potential and thus, provided the changes 
are made sufficiently slowly, will not change the universal aspects of the coarsening provided all lengths are scaled 
appropriately. 

The most obvious difficulty for any experiment is one that is ubiquitous in random systems: how does one obtain 
a reasonable range of length scales if the correlations are growing only logarithmically in time? The situation here is 
not quite as bad as it might seem as on macroscopic time scales £ <~ In 2 t is not all that short: with a microscopic 
time, ro, of order picoseconds, a minute corresponds to coarsening by a factor of a thousand in length scale. But for 
one-dimensional random-field systems, one can do much better: the correlation length is of order 



with ao a microscopic scale. Thus the correlation length can be controlled by a combination of increasing time, 
increasing T and decreasing h. To get a wide range of £(t), it is beneficial to start with h large enough that 
the shortest time scales which can be measured still correspond to microscopic lengths and then decrease the field 
gradually. 

While some of the non-equilibrium properties studied here should be amenable to conventional experimental probes 
in magnetic systems, some of the most interesting predictions will probably not be: how the specific set of domain 
walls depends on equilibration time and how it varies — or does not vary - from run-to-run on the same sample. For 
these one needs either microscopic imaging probes — perhaps with magnetic atomic force or tunnelling microscopy — 
or scattering data with enough resolution that speckle patterns from a finite spot size can be measured. As mentioned 
earlier, the spatial Fourier transform of the two-point two-time correlation calculated here is, for spin glasses, related 
to the correlations between the speckle patterns at the two different times. While measuring this with magnetic x-ray 
or neutron scattering may not be possible, it should be feasible via light scattering on systems in the same universality 
class in which the length scales are sufficiently long. 

One system on which light scattering measurements might be possible is a nematic liquid crystal in a long thin 
tubes with a square crossection and heterogeneities on the surfaces which would couple randomly to the two possible 
symmetry related orientations of the nematic director. Whether this or other systems can be formed in a regime in 
which the dynamics are sufficiently fast to allow a wide range of length scales to equilibrate is a question whose answer 
must rely on a quantitative analysis of the physics on the scale of the domain wall structures that should occur. 



Many of the qualitative features of the coarsening process in RFIM chains are also expected to obtain in higher 
dimensional random systems. These are particularly interesting in systems in which, in contrast to the one-dimensional 
models, true long-range order should exist in equilibrium. Nevertheless, the characteristics of the coarsening process 
in the one-dimensional models with weak randomness should have much in common with such systems as long as 
consideration is restricted, as we have primarily done, to time-scales such that the full equilibrium correlation length 
cannot be attained. 

Features that are common to many higher dimensional random systems — random exchange ferrommagnets, random 
field magnets and, although more controversially, spin glasses — include: growth of some kind of domain structure by 




(246) 



B. Higher-dimensional systems 
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thermal activation over barriers that grow with length scale and are broadly distributed; and the existence of aging 
and other history dependent phenomena. From a renormalization group point of view, these features are general 
characteristics of systems controlled by random zero-temperature fixed points [5,41]. The notions of local equilibrium 
within constraints caused by large barriers and of domination of the dynamics at any given long time scale by rare 
regions of the sample in which the rate for going over barriers is of order the frequency being probed, are both 
important. In general, this will, as in the one-dimensional models, result in three different regimes for two-time 
properties depending on whether the difference between the times is much smaller, of order of, or much larger than, 
the earlier time. 

One of the most intriguing questions concerns the pseudo-determinism found in the one-dimensional random models. 
Will this exist in some higher dimensional systems in which the domain walls arc lines or surfaces whose evolution in 
time will involve topological changes, rather than points which can simply move or annihilate? I.e., will domain walls 
tend to lie, at long times, in similar sample-and-time-specific positions which are only weakly dependent on the initial 
conditions? Or will they evolve in different runs in very distinct ways- perhaps by retaining much more memory of 
the initial conditions? 

Two examples of a single random walker diffusing in random potentials illustrate some of the difficulties of drawing 
any definitive conclusions. Consider, first, a random potential which is independently random at each site with a long 
power law tail ot the distribution of the depth of the potential wells. The deep wells will give rise to subdiffusive 
behavior, indeed, simple considerations of the time to find a deep trap and the time to escape it imply that the typical 
distance the walker will diffuse in time t only grows as a power of Int. But the statistics of the set of sites the walker 
visits by the time it has made a given number of jumps from one site to another will be identical to that of a normal 
free random walker; it is just the time spent on each site that causes the slow diffusion. Since a pair of random walkers 
in dimensions higher than four have a non-zero probability of never visiting any sites in common, it is clear that the 
long time behavior is not deterministic in high dimensions: two different runs starting from the same point will have 
a probablity that vanishes at long times of the two walkers being at the same site. In two dimensions, in contrast, the 
pair of walkers will be very likely to be trapped at the same site at long times. [The three dimensional case is more 
subtle but will be more like the high dimensional that the low dimensional case.] 

A second type of random potential yields different behavior: a gaussian random potential with mean zero and 
[V(x) — V(y)] 2 <~ |x — yl 2 '' with £ > 0. A pair of random walkers in such a potential will tend to get trapped in the 
same deep valley at long times. This can be argued by considering the borders of valleys which are surrounded by 
barriers of at least a height T. If these typically do not have a lot of fine-scale structure, then they should act as 
effective traps for all walkers in the vicinity. All the walkers will then tend to eventualy leave such a valley over the 
same barrier and hence end up in the same next-larger-scale valley as well. The behavior would then be asymptotically 
deterministic as in the one-dimensional Sinai model. 

The crucial feature that distinguishes these two cases, seems to be the smoothness of the potential: in the first 
example, the deep potential wells come from very short wavelength fluctuations and there is thus no useful notion 
of a coarse-grained potential. But in the second example, the deep potential wells come, as in one-dimension, from 
long-wavelength slowly varying components; thus coarse-graining the potential to yield approximate dynamics which 
are asymptotically exact should be possible. We should emphasize, however, that it is by no means established even 
for this simple model of a random walk in a smoothly varying gaussian random potential that the conclusion argued 
above — that the dynamics will be asymptotically deterministic — will be valid in high dimensions. 

What happens in three dimensional random magnetic systems in which there are truly many degrees of freedom, 
we must leave as an important open question. 
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APPENDIX A: CONVERGENCE TOWARDS THE FULL STATE 

In this Appendix we analyze the rapid convergence towards the full state that was discussed in Section (IV C). 

We consider, for simplicity, random initial configurations of the spins, in which there is an independent probability 
w that each intersite position is occupied by a domain wall, and a probability 1 — w that it is not. Such random 
initial configurations describe, for example, initial equilibrium before a quench from a temperature T which is 
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high enough that the random fiels are negligible and the system behaves like a pure chain. This corresponds to 
T = l//?o {h n } and the choice w = 1 ^. e -^ 2j ■ In particular, initial conditions where all spins are independent 

and take value ±1 with probability 1/2 corresponds to infinite temperature, /3 = 0, and w = 1/2. 

For these type of initial configurations, the probability that there exist n = 0, • • • I domain walls among I consecutive 
lattice spacings is simply given by the binomial distribution 

^I0 = ^T-"(1 -»)'-»• (AD 

In the renormalized landscape, the length I of descending (resp. ascending) bonds at scale T is distributed with 
-Pp"(Z) (resp. Pf(l)), whose Laplace transforms are obtained after integration over £ in (37). The probability that 
there had been initially no domain walls in the interval occupied by a descending (resp. ascending) bond of the 
renormalized landscape plays an important role; it is 



zero 



(T) = J2 P r(l)RW) = E P r (0(1 - w) 1 (A2) 



1=1 1=1 
The probability that there had been an odd number n is 



rU?) = E P r (0 E *M0 = E P r(0 1 ~ (1 ~ 2W) (A3) 

/— 1 n odd 1=1 

and the probability that there had been a non-vanishing even number is 

rfven (r) = 1 - r± d (r) - r± ro (r) (A4) 

To characterize the statistical properties of the spin configurations on the renormalized landscape, we focus on one 
renormalized bond at scale T — call it "2" — and its local environment which determines whether or not the maximum 
and minimum at the ends of bond 2 are occupied by domain walls. We introduce the probabilities V^ esc (e2] {ei, £3}) 
(resp. Vp lsc (e2; {ei, £3})), where = ± , that a descending (resp. ascending) bond of the renormalized landscape at 
scale r is in phase ti = ±, with its left neighboring bond in phase t\ and its right neighboring bond in phase 63. The 
normalisation conditions are 

E E E V r d ^(e 2 ;{e 1 ,e 3 }) = l (A5) 

ei=± C2=± f-3=i 

E E E VT c (e 2 ;{ ei ,e 3 }) = l. (A6) 

ei=± C2=± C3=± 

If e 2 = — on a descending bond, then it is not possible to have ei = + since this would correspond to a domain 
wall of type A at a maxima, and similarly it is not possible to have e 3 = + since this would correspond to a domain 
wall of type B at a minima and thus we have immediately that 



Similarly we have 



V r desc H {+-}) = y r desc (-; {-+}) - V r desc (~; {++}) = (A7) 



V r asc (+; {+-}) = V r asc (-; {-+}) = V r asc (-; { — }) = (A8) 



If, however, e 2 = + on a descending bond (resp. f2 = -ona ascending bond), then all four possible neighborhoods of 
this bond can occur. Since it is a bit lengthy to give the full enumeration of all possible cases with their corresponding 
probabilities, we give here only the final results 



V desc(_. 1 }) = 'zerJT) (AQ) 

^r sc (+; {++}) = r * er ° (r) (A10) 



and 
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V 4esc {+ . { __ }) = (1 _ r - ro(r)) ^ _ rJer^T) r+ ro (r)r~ ro (r) ^ + (r~ r .(r)) 2 ^ (r) (An) 

^ 8C (+;{+-}) = U r desc (+;{-+}) = - (rzer ° 2 (r))2 (r+ ro (r) + r+ eTt (r)) (A12) 

U r desc (+; {++}) = (r -° 2 (r))2 (r+ ro (r) + r+ e „(r)) (A13) 

and similarly for ascending bonds by exchanging ± — > =p. 

The important property of these probabilities is that apart from V^ esc (+;{ }) and V^ sc (—; {++}), which cor- 
respond to locally full configuration of the domain walls, all the allowed V have rf ero (T) as a factor, i.e. to have a 
bond (2) that does not have domain walls at both its extremities requires that at least one of the three bonds (1,2,3) 
had exactly zero domain walls in the initial configuration. 

Using the fixed point solution (37), we find that 



zero 



(I> dlP±{l){l - w) 1 =P?(q = In (j^h;)) (A14) 



and thus this and the probabilities of non-full bonds converges exponentially in T to 0: 

rt ero (T)~e-^V^^) (A15) 

withp=^ln( T ^) 

Thus the system converges towards the "full" state of the renormalizcd landscape exponentially in T, with a non 
universal coefficient depending on the initial concentration of domain walls through the parameter w, and on the 
strength of the disorder through g. For the symmetric case of no applied field (5 — 0), the fraction of the extrema of 
the renormalized landscape at time t that are not occupied by a domain wall goes to zero as a power of time: 

Prob [missing domain wall] <~ — (A16) 

with 



(A17) 

Note that in the initial state for the landscape, and in the symmetric case H = 0, h n is positive or negative with 
probability 1/2, and thus after grouping the consecutive h n of the same sign to give the initial zig-zag landscape, we 
have for initial distribution of length 

Po(l) = ^ for J = l,2,---oo (A18) 



This corresponds to 



,(r = 0) = l^ (A19) 



As an example, for w = \, we have r zero (T = 0) = |, and Vfl s °(—; { }) = |. 

APPENDIX B: AUTOCORRELATIONS IN THE RFIM WITH AN APPLIED FIELD 

We have now to solve the RG equations (94) given in the text with the initial condition 

p++,(o= r^/^M (B i) 

JO IT' 

P-+,(C)=0 (B2) 

Since for large T, P ± have reach their fixed point values (37), we obtain ( using the simplified notations Up = 
4(p = Q) = U±(p = 0) ) 
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r°° + + 

(d r d c )P±+(() = -2ufP r ± +(C) + 2u*u± / dC' e -"r(C-C') P ±+ ( + (4) 2 Ce-"rCpT+( ) (B3) 

Jo 

together with the initial conditions 

P++,(C) = =Le"^ (-a p C/ r +(p) + CU+(p)d P u+(p)) Ip=o (B4) 

P r -+,(C)=0 (B5) 
The solutions are thus of the following form 

P±+ (C) = + (B±+) (B6) 

where the coefficients satisfy 

p±+ = a r A±+ (B7) 

<9;U±+ = u+Mp + Ap + ) (B8) 

Introducing the sum £p r , = Ap p/ + Ap p, and the difference Ap r , = p, — App", , we obtain the decoupled equations 

^A+p, = (B9) 

^S+p, = 2M+MpE+ r , (BIO) 

together with the boundary conditions 

A+ r , = S+ r , = = ^ (V coth 7 ' - 1) (Bll) 



drA+ r ,|r=r" - d r Z+ rl \ r= r> = B++ r , = — coth 7 ' - , 1 - ) (B12) 



^ r '-^( COth7 '-SnnV 
where 7 ' = <5r" as usual. In terms of the function 

p(j) = 7Coth 7 - 1 (B13) 

and its derivative with respect to 7: 

M{i) = coth 7 ^— ( B14 ) 

sinh 7 

we finally obtain 

A rP = ^2 ± 7M(Y) ± P(Y) T 7'M(Y)) (B15) 
B±+ =1(M(7)±M( 7 ')) (B16) 

In the same way we obtain the solutions for the minus initial condition 

P±p-,(C) = =e~ C4 4 + <B^) (B17) 

with ^4p p, = Ajfp, and Pp p, = Pp p, . We may check the normalizations 

l^dC^P+Uo + Pf^o) = i(i± M{i)) (bis) 



The autocorrelation is 



< ^ (f )5i (t) > = l°° d ( ( p r.r (0 + ^ (0 - ^rtr< (0 - ^ (C)) (B19) 

^ 1 (2 - Ap" r ,) + (-L + ^) - Pp- ) ) (B20) 

= (coth 7 )M( 7 ') + -rL- ( 7 M( 7 ') + P (i) - l'M{i)) (B21) 
sinn 7 



Note that < S^t') > = M{i) and < S l {t) > = M(j) = coth 7 - 
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APPENDIX C: TWO-POINT TWO-TIME CORRELATIONS {S (T)Sx(T))(S»(T')S x (T')} 



1. Definitions 



Two points at and x > are given. To compute the two spin two time correlation we introduce the following 
quantities. Let np" r 7(Co,2o; Ci,Zi, C,2n, hn', x), n = 0,1,2.., (and respectively Qp" r 7(Co,Zo; (i,h, ■■■', C2n, hn', x)) be 
the probability that at V , x and are on bonds of same orientation (respectively different orientation) AND that at 
r, x and are on bonds of same orientation, with a configuration (Co, W, Ci, h, ■ ■■', C,2n, hn) (sec figure). One defines 
similarly f^r?r' 1 ' + (Cn, ^o; Ci, C271+1, hn+i, x), n = 0, 1, 2.., (respectively fip™^ 1 ' - (Co, ^o; Ci^i, Can+i, hn+i', x)) be 
the probability that at V , x and are on bonds of same orientation (respectively different orientation) AND that at 
T, x and are on bonds of different orientation, with a configuration (Co, ^o! Ci, h, ■ ■■', C2n+i, ^n+i)- 

Initial condition at T = T' (for all n > 0): 

^r^rKCo, h', Ci, h, (2n, hn', x ) (CI) 

2n 



= P v , (Co, h)Pv (Ci.M--- P T> (C2n, hn)W T> (l , hn, ^ h - x) (C2) 

i=0 

fi r ?p/' (Co, W, Cl, h, C2n+1, hn+1', x) (C3) 

2n+l 

= Pv(Co,h)PT'(Cl,h) ■ ■ .Pr>((2n+l,hn+l)Wr>(l ,hn+l, Yl h - x) (C4) 



with the notation 



w r (h,h, £) = 4 t /' 2 d y^( L - »i - 2/2) (C5) 

1 Jo Jo 

2 

= — [mm(Zi, L) — mas(0, i — h)\ [min(li, L) — max(0, L — h)\ (C6) 
We have of course Qp?'p, = fip?^, 1 ^ = 0. 

In the end we are interested in the probabilities Pp'p, (x) that (So(t')S x (t')} = e' and (So(t)S x (t)) = e. We have the 
normalisation 

P+P+ (x) + Pff (x)+ p-p+ (x)+ Pffi (x) = 1 (C7) 

The correlation function is 



(S (t)S x (t))(S (t>)S x (t>)) = P r + r t(x) + P r : r 7(x) - Pffi{x) - Pffi(x) (C8) 
In terms of the functions Q we have 

Opp7(Co, h', Cl) h, C,2n, hn', x) (C9) 

^r™f" (Co, Cl) 'l> C271, hn', x) (CIO) 

Q r ™+ ' + (C , h; Ci, h, C,2n+i,h n +i',x) (Cll) 

■.+ I,i2n + 1 

^p 1 ' (Co, ^o! Ci, ^i) •••) C2n+i) ^2n+i; a;) (C12) 



Pp+ r t(x) 


00 „ 






Pp+ r 7(x) 


00 „ 




„=0 J Co,2o, 


P r - r t(x) 


00 „ 
n ^ ^Co,/o, 


Pp r , (x) 


CO „ 

,i=0 J| > ''°> 
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2. RG equations 

The RG equations for Q™'^, for e' = ± and m = 0, 1, 2, . . . read : 

I _ X] ^ fc ) ^™r' (Co, ^o; Cii'i; • • ■ Cm: l m ;x) = 

\ k=o ) 

-2P r (C = 0)fi™r e ,'(Co, Z ; G, ■■■Cm, l m ; x) (C13) 

TO „ 

+ X] / f^p", 2,6 (Co,Jo; • ■ • ; Cfe-i, ife-i; ^, Z; 0, 1'; Cfe - 2, Z"; ■ • ■ ; CmJm,;x) (C14) 

k=0 z - l+l'+l"=lk 

+ [ P T (Co-z,l)^+y(o,l' ;z,l' 1 ;Ci,h;...;C m ,l m ;x) (C15) 

Jz,l+l' +l' 1 =l 

+ Opp, 1 ' 6 (Co,io;---;Cm-i,im-i;^,C;°'C+i; a; )- F i , (Cm-2,Z) (C16) 

+ [ p r (z,i)P r (o,i')n^' l (c -z,i' ;Cuh;...;Cm,i m ;x) (Ci7) 

+ / ^r.r' (Co ; Zo! • ■ • ; Cm— ij Z m — ij z, Z m ; x)P r (0,Z')Pr(C P -^0 

Jz,l>+l'+l=lm 



In particular, for m = we have 



(C18) 



(9r - d Co ) Op p, (Co, la;x) = 

-2Pr(< = 0)Qp' e r ,(Co,Zo;z) (C19) 



+ / fi r ' e p,(z,Z;0,Z';Co-^,Z";x) (C20) 

Jz,l+l'+l"=l 

+ [ Pr(&-z,l)i%i,(0,k,z,l' 1 ;x) (C21) 

+ / ^ r ; e r,(z,Z";0,Z';x)P r (Co-2,Z) (C22) 

+ [ Pr(z,l)Pr(0,l')n r ^,{<:o-z,l' Q ;x) (C23) 

+ / r! r ' e r,(z,Z o ;x)Pr(0,Z')Pr(Co-^Z) (C24) 

J z,l' +l'+l=l„ 

+ [ Q°i,(0,l' ;x)Pr(z,l')MCo-z,l) (C25) 

Jz,l' +l'+l=l 

3. Form of solutions 

For n > 1 and e' = ±, we set 

2n 

^r"f' (Co, Zo; ■ • ■ ; (2™, '2™; a;) = Pr.r'(Co, Zo; C,2mhn\ k — x)Pr(Ci,h) . . . Pri&n-iihn-i) (C26) 

i=0 

2n+l 

fip™^ 1 ' 6 (Co, Zq; ■ • ■ ; C2n+1, Z2n+l! X) = E T p, (Co, Zq; C2n+1, Z2n+l! ^ Z, — x)Pr((l, M • • • Pr(C2ri, Z2n) (C27) 



i=0 

and also 



^r',r'(Co,Zo;Ci^i;a;) = ^r,r'(Co,Zo; Ci^i^o + Zi - a;) (C28) 
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where the Pr(C>0 satisfy the bond equation, and where the E satisfy the RG equations 

(dr - % - %)i^;(Ci, h; ( 2 ,h;L) = -2P r (( = 0)^'(Ci, h; C2, h; L) (C29) 

+P r (0, .) * h Pr{., ■) *a,h E±f',(., .; C2J2; L) (C30) 

+P r (0, .) * h Pr(., •) * C 2,h E±£',((uh; •, •; £) (C31) 

+ / £; r F r e ;(0,r i ;C2,Z 2 ;i-0^r(.,0*C 1 ^r(.,n (C32) 

+ [ ^,'(Ci,/i;0,^;L-0i¥(.,0*c a ^(-.O ( C33 ) 

+ / £ r t ^'(.,Z / 1 ;C2^ 2 ;i-/-0 * Cl P r (.,0^r(0,O (C34) 

+ [ £; r t ! r e ;(Ci,Zi;.^ 2 ;L-?-n*c 2 ^r(.,0^r(0,n (C35) 

J2+i'+i 2 =Z 2 

with the initial conditions 

£ r ',f r ,(Ci, /1; C2, fe; L) = 5 e , e ,Pr>{Ci,h)Pr>{C2,l2)W r >{li,h, L) (C36) 

The RG equation for O r 'p, becomes 

(dr - d Co ) n°£, (Co, W, x) = -2P r (C = o)n°£< (Co, l \ x) (C37) 

+ / £+; r <OM;Co-^'Vo-^r(0,0 (C38) 

J z.l+l'+l"=la 



with the initial condition 



(C39) 



+ / P r (.,0*Co4',r<(0,4-,Zi;z) 

+ / n^r,(.,40,Zi;x)* Co iV(.,0 (C40) 

+2P r (0, .) * ;o P r (., .) * Co , io .;*) (C41) 

+P r (., •) * C o,io *H> •) * io n°£,(0, .; x) (C42) 

^ACo,lo\x) = 8 e > t+1 P r >(Co,lo)W r >(lo,lo,lo-x) (C43) 
Using the form of the solutions, we thus obtain 

P+; r +(x)= [ Q°'+(Co,Zo;z) (C44) 

00 . . 2n 

+ ]T / / ^+(Co,/o;C2r l ,^;i)<5(L-(^? l -x))P r ai)...Pra 2 n-i) (C45) 

Prf(x) = f n r '-,(Co,W,x) (C46) 

00 „ . 2 n 

+ ]T / / E+^(( ,l -,<; 2n ,l 2n -,L)6(L-(£2k-x))Pr(h)...Pr(l2n-i) (C47) 

Pf${x) = (C48) 

00 . „ 2n+l 

]T / / Pp>+(Co,?o;C 2 n+i,Wi;i)^-(E ^-^r(/i)...Pr(/ 2 „) (C49) 

„ = ^Co>C2n + l J L,l ,h,l2n+1 i—Q 

Pff(x) = (C50) 

00 „ „ 2n+l 

53/ / £p>" (Co, *o; C 2 n+i,W; (E ^-^)) p rGi)...Pr(M (C51) 

n=0 ^Co,C2n + i J L,l ,ll,l2n+1 i = Q 
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4. Laplace transforms 

It is convenient to introduce the Laplace transforms 

/>oo />oo />oo 

£r'r<(Ci,Pi;C 2 ,p 2 ;p)= / d*i / di 2 / ^e- piZl - p2Z2 - pi i?r.r'(Ci^i; C2, h\ L) 
Jo Jo Jo 

Using the fixed point solution 

P r ((,p) = U r (p)e-^(p) 

the RG equations for E become 



(dr - % - %)^'(Ci,Pi;C2,p 2 ;p) = -2^ r (0)^;(Ci,Pi;C2,P2;p) 

+UKpi) dze-^- z ^^E^:(z, Pl ;C 2 ,p 2 ; P ) 
Jo 

+U*(p 2 ) [ 2 dze-«>-'*rWE±g(C uPl -,z,p2-,p) 
Jo 

e -Ci«r(pi+p) _ g-Ciur(pi) , 

+U r (pi ) U T (pi + p) — , — — ; 3 ^ (0, pi ; C2 , P2 ; p) 



+£/r(p2)?/r(p2 +p) 



wr(pi) - ur(pi +p) 

e -C2«r(p2+p) _ e -C2"r(P2) 



^r,r'(Ci,Pi;0,P2;p) 



ur(p 2 ) - u r (p2 +p) 

+UKP1 + P) C 1 dze-^-^^+riE±'*',(z, Pl ;( 2 , P2 ;p) 
Jo 

Jo 

with the initial conditions 

2 [j7r'(Pi)e~ Cltlr ' (pi) - t/r-(pi +p)e" ClMr ' (pi+p) 
t/ r '(P2)e~ C2Ur ' (p2) - U r >(p2+p)e- C2Ur ' ip2+p) 



Er', e r ,(Ci,Pi;C2,p 2 ;p) = $e 



Also defining 



^r"r°' £ (Co,Po;p) 



rOO rOO 

Jo Jo 



the RG equation becomes 



(dr - d Co ) Op'p,(Co,Po;p) = -2u r ^r ; r/(Co,Po;p) 
+ / ^r r,(2,p +p;Co-2,_po+p;-p)^r(0,po+p) 

+(7 r (p ) [° dze-«°- z ^^n r i,(0,p o -,z, Po ;p) 
Jo 

+Ur( P o) [ ° dze-^- z ^ ( ^n r i,(z 7 p o ;0,po;p) 
Jo 

+2f/2( po ) /" ° dze-^- z ^^n r i,(z,p ;p) 
Jo 

+U^po)Coe- CoUripo) n r ^,(0,po-,p) 



(C52) 

(C53) 

(C54) 
(C55) 

(C56) 

(C57) 

(C58) 

(C59) 

(C60) 

(C61) 
(C62) 

(C63) 

(C64) 
(C65) 

(C66) 

(C67) 

(C68) 
(C69) 



with the initial condition 
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T ,2 p 



-pd P0 (t/r'(po)e- Co " r ' (po) )] 



In Laplace 



e 4,{q) = / dxe- qx P^i,{x) 
Jo 



the relations (C51) become 

p+i+(g)=/ n°-+,(Co,o ;g )+ - Pr ^) / s++(Ci,g;C 2 ,g;-g) 

p&{q) = J ftf>(Co,0;<?)+ 1 ^ 1 ^ S+-,(Ci,g;C 2 ,g;-g) 

p r\r!"(9) = i / £ r,r'(Ci> 9; C2, g; -g) 

1 - Jd,C2 



p r,r' (g) = ! _ p2( g ) ^ ? £ r>r ,(Ci,g;C2,g;- 



5) 



(C70) 
(C71) 

(C72) 

(C73) 
(C74) 
(C75) 
(C76) 



Thus for the functions S, we need to solve only the case pi = p 2 = q = —p, i. e. with the notation Ur(p = 0) 
u r (p = 0) = u r 



{dr-d (l - d C2 )E^ f,(( l7 q-X 2 ,q; -q) = -2u r E^, (Ci, q; C2, g; -q) 

+U*(q) dze-&-'*rtoEpg(z, g; ( 2 ,q; -q) 
Jo 

+U*(q) f * dze-^- z >^E±i{Q liq] z, q --q) 
Jo 



+U r (q)u r - 



-Ci«r _ P -Ciur(q) 



-E^',{0,q;( 2 ,q;-q) 



u r (q) - Mr r ' r ' 

e -C2«r _ e -C2"r(g) , 

+Ur(l)ur ^773 — — (&, «; 0, g; -9) 



ztr(g) - w-r 

. dze-^-^Ep^(z iq -X 2 , q --q) 
Jo 

+4 [ 2 dze-^- z ^E^:(( u q;z,q;-q) 
Jo 

with the initial conditions 

2 



£ T'T'(Cl,<7;C2,<7;-<7) = <Se,e 



p2 g 2 



[7 r '(g)e- Cl "r'(9) - Uvie -Ciu T , 



For n = 0we need to solve only the cases po = p\ = and p = q, i.e. the equation 

(dr - d (o ) (Co, 0; g) = -2« r ^ (Co, 0; g) 



+C/ r (g) / E^,(z,q;Co- z,q;-q) 

J Z 

Jo 

+u r [ dze- (c °- z)ur f7 r 'p,(z,0;0,0;g) 
Jo 



(C77) 
(C78) 

(C79) 

(C80) 

(C81) 

(C82) 

(C83) 

(C84) 
(C85) 

(C86) 
(C87) 

(C88) 

(C89) 
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+2u 2 T [ ° dze-^°-^ ur n%(z,0;q) (C90) 
Jo 

+u 2 r ( e- CoUr n° r ^,(0,0;q) (C91) 

with the initial conditions 

n^' r ,(Co,0;g) = 5^ +1T ^lU r ,(q)e-^>(i) _ £/ r , e -<o«r< (C92) 

-pd P0 (U r > ( P o)e'^' ^) | Po=0 ] (C93) 



5. Solution for the functions E 



Using the symmetry in (£1, £2), we set 

^r'(Ci,9;C2,«;-9) = ^(4;r-(?)e- a,,r(,) e- C2UrW + B^^je"^^-^^) 

It is in fact more convenient to use the following combinations 

r 2 

FtHv) = Y^r'r'(Ci = 0, q; C2 - 0, q; -q) 
= A^,(q) + 2B^,(q)+D^,(q) 

r 2 

d( 2 Erf r , (Ci = 0, q; C2, q; -q) 



Jr'r,(q) 



K?Uq) 



2 jo 

A rr (g) 
ur(q) 

r 2 
2 



/•CO /"C 

/ dc w 

Jo Jo 



+ B^ rl (q)[^ T - + — 



y r,r' 
Mr 



rfC2-E r ' e r'(Ci,9;C2,9; -<?) 



A r'r'(«) , „ B r.r'(9) , ^r.r'te) 



Up(q) u r (q)u r 



They satisfy the system 

0r^£'(g) - -2(u r + w r (9)Kr,r'(9) + 2« r «r(«)^/'(g) 



^ «r(«) > ) ^'(9)-^'(9)+«r«r(«)^'(9)+ ' ' ^ 



and the initial conditions 



«r(g) 
'ur(q) 



ur(q) 



2u T ) K^,{q) - 2J^, (q) + 2^f J+f, (q) 



T±,e' 



Ur(q) 



r/ r ,(q) = 6 e , e ^[Ur'(q)-ur'} 2 



Ur>(g) 
ur'(q) 



K^ T ,{q) = 5, 



q 2 



Ur>(q) 
«r'(,) 



1 2 



It is convenient to introduce the following sums and differences 



(C94) 
(C95) 

(C96) 
(C97) 
(C98) 

(C99) 
(C100) 

(C101) 

(C102) 
(C103) 

(C104) 

(C105) 
(C106) 

(C107) 
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4\r> (<?) = 

4'Uq)- 

K^,(q) 



Ip^,(q) - ip^,{q) 

= J+4,(q)- Jr4,(q) 
= Kp$,( q )+K-^,(q) 
= KpUq) - K-'Uq) 



Then for e' fixed, the three functions I S e , J S ' e and K S ' e satisfy the system 
9r/p'p,(g) = -2(u r + ur{q))4'^,{q) + 2uru r (q)4'r'{q) 



l r,r' 

T S,e' 



d r J r ' r ,(q) 



' r( '" -u r (q)) 4i(q)-4i(q)+u r u r (q)K^(q)+ l ^4i(q) 



ur(q) 



drK?i(q) = (2^1 + 2u r ) K^,(q) - 2J^(g) + 2^J& 



(9) 



with the initial conditions 



/ r ; ,r'(9) = ^[ f/ r'(«)-^] 2 



4x> (1) 



1 



[U r >(q) - u r , 

Ur>(g) 
u r >( q ) 



Ur-{q) 
u T '{q) 



- 1 



1 2 



and thus the solution is simply 



4:ki) = i[Mi)-UTY 



J r,r'(l) = ^[^r(g)-«r] 



<'r'(?) = 4 



u r (g) 



- 1 



For e' fixed, the three functions I D ' e , J D ' e and if ^ satisfy the system 
d r Irx' (l) = -2(?i r + u r (q))I^, (q) + 2u T u T {q)J°£, (q) 



u r(q)) Jr.r> (<?) - Ir.r' (l) + «r«r(?)^rr' (?) _ ~^T^ 7 rr' (<?) 



«r(g) 



Sr#££ (q) = 2 



"r : '/; 



(C108) 

(C109) 

(C110) 

(Clll) 

(C112) 

(C113) 
(C114) 

(C115) 
(C116) 

(C117) 



(C118) 
(C119) 

(C120) 

(C121) 
(C122) 

(C123) 

(C124) 
(C125) 

(C126) 



with the initial conditions 



I?4,(q) = (-iy'±lU r ,(q)-u r ,} 2 
J^,( q ) = (-iy'l[U r ,(q)-u ri } 



Ur'(q) 
u T >{q) 



Ur>(g) 

Ur>(q) 



- 1 



(C127) 
(C128) 
(C129) 



Three linearly independent solutions of the system read 
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I» 1 (q) = {U r (q)+u r y 
J° 1 (q) = {Ur(q)+u r ] 



Ur(q) 



K?\q) = 



Ur(q) 



«r(,) 



(C130) 
(C131) 

(C132) 



I? 2 {q) = 2u r 



Urjq) + ur 

U r (q) + u r (q) 



JF 2 (q) = 



(2u r + U r (q)) 



Ur(g) 
«r(<?) 



K?\q)=2. 



Ur(g) 



Ur(q) + ur(q) 
1 



wr(g) 



(C133) 
(C134) 
(C135) 



r C3 



(q) 



Jr D3 (q) = 



KP 3 (q) = 



Uv{q)[2u 2 T + (2u r + Efr(g))("r(g) + Ur(g))] 

(ur(q) + U T (q))(u r (q) 2 -U r (q)*) 
U r (q)[2u r + u r (q) + U r (q)} 
«r(9)( Wr (g) 2 -C/r(g) 2 )) 

2U r (q) 



U T {qf{u T {q)-U T {q)) 
It is useful to consider the matrix formed by these solutions 



(C136) 
(C137) 
(C138) 



Nr 



ifHq) 

-Dl 



(q) 



IfHq) 

jD2 
J r 

D2( 



(q) 



If 3 (q) 
JfHq) 



V K»\q) K™(q) K^(q) 



I 



whose determinant gives the Wronskian Wr of the three linear independent solutions: 

^ r = det[7V r ] = -^fM 

The solution satisfying the initial conditions we are interested in will be obtained as a linear combination 



(C139) 



rD,e', 

'r,r' 



V K^',{q) 



\ 



J 



/ iP l {q) 

jPKq) 



i=i 



V 



Kp\q) 



) 



\¥(q) 



\ 
J 



where the coefficients X r e , (q) are determined by the initial conditions 



\ T i{q) 



KP4'(q) 



The inverse of the matrix Nr is 

/ 2u T (q) 



I 



(C140) 



N r 



V 



Ur(qy 2 [u r (q)-Ur(q)] 

(ur(q)+Ur(q)) 
Ur(q) 2 

lur(q)+Ur(g)] 2 
Ur(q) 2 



2u r (q)(2u r ur(q)+Ur(q)lur(q) + Ur(q)]) 

u v (q)(2u T +U v (q)) 
Ur(q) 2 

2uy (q) [ut + Ut (<?)] [ur (q)+Uv(q)} 
Ur(q) 



u r (q) 2 (2u 2 ,u r (q)+U r (q)[2u r + Ur(q)][u r (q)+Ur(q)]) 
Ur{q) 2 [ur(q)+Ur(q)][ur(qy-U r (q) 2 } 

u r u r (q) 2 [u r + U r (q)] 
Ur(q) 2 [ur(q) + U r (q)] 

ur(g) 2 [ur+Ur(q)] 2 
WW 



and finally we get the coefficients 
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,u' ( , U e' (-8u*,u r ,(g) + Mg) + ^ r -(g)][-5^(g) 2 + E#(g) +M3Mg) - ^(g))]) 
r ' (<Z) " ^ (ur , ( , ) _ [/r((g))(ur/(g) + [/r , ((z))2) (C142) 



A^( g ) = -(-l)'' 4(ur '-r (g))2 (C144) 
? 

In the end we are interested in 

P&(q) = J ^(Co,0;«) + f^ i ^2^) ^(9) (C145) 

^r> (?) = jf ^r>(Co,0; 9 ) + r2 ^ t ^p2( g) g r,r" (g) (C146) 

^(g) = ^ i_^ (g) g r,^(g) ( C14? ) 

p r~> (?) = ^ i_p2 (g) ^F,b"(g) ( C148 ) 

Since we have the relations 

K^+(q) = K^-(q) (C149) 

K°$(q) = -K°£(q) (C150) 

(C151) 

we can express the four functions K^ e r ,(q) in terms of K r '^, (q) and K®\5(q) only as 

Kr'Uq) = \ «£(?) + (-l) £ (-l) e '< r +(g)) (C152) 

Since we have the following constraints 

P+i+ (?) + P r + r 7 (?) = P r + (?) (C153) 

P r "i+ (q) + Pffi (q) = Pf (q) (C154) 

P r + r t (?) + P r ~>t (?) = Pft (?) (C155) 

P+ r 7 (?) + Pff (?) = P r " (?) (C156) 

where 

*<«> ^ f^2 £gfj£ (OK) 

we do not have to compute separately the functions Qp p,(Co, 0; ?). 
Indeed, we get 

P^(q) = J f2 r '^(Co,0;?) + ^ O r ; r ,(Co,0;?) + A _j^|_ ^ ; + (g ) (C159) 

^r-(9) = ^ i_p2 (g) ^r>(g) (C160) 

= / ^'(Co,0;?) + ^ ^((1 + Pr(?))<'r",(?)-(1-Pr(?))^ r t(?)) (C161) 

Pf>{q) = J ftr>(Co,0;?) + !_p2 (g) ((1 + Pr(?))<'^(g) + (1-Pr(?))^rt(?)) (C162) 
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The second equation is satisfied since K^,(q) — \ [1 — Pr(q)] z . The three other are compatible and give 



/ 

L 



« fJ ,( (o ,0; 9 ) = i - ^ - J^d - i*( 4 )) + j^Jf^ta) (C163) 

«f - r ,( C „,0;,) = jA, - J^d - «.(,), - i; T^ASth) (C164) 



6. Final result 

The Laplace transform of the correlation function can now be obtained as 

/ dxe-« x (S (t)S x (t))(S (t>)S x (t>)) = P+i+(g) + Pff(q) - Pfj+(q) - P+f(q) (C165) 
Jo 

K°;+{q) (C166) 



1 4 , 4 1 T ^D, 



q T' 2 q 2 1 + Pr,(q) T 2 1 - Pg(q) r > r 
and thus the final explicit expression is 



dxe-" x {S (t)S x (t)){S (t')S x (t')) = - q - p^tanh 2 (^j^j (C167) 
+ ±cotanh 2 f X r t(q) + 4 ^|^ A 2 t( 9 ) + 4 . ^ ^(g) (C168) 



V 2 ; r w r 2 r ' w r^ sinh 2^ 



where 



X r t(q) = 2r^ g 3 si nh 2 (r'Vg) [8 + 3T ' 2q ~ 16cOSh(r '^ } + 8 r 'V^sinh(r'^) (C169) 

+ (8 + 5r' 2 g) cosh(2r'V?) - IW'Vq sinh(2r'Vg)] (C170) 

X T t(q) = ~^{fi- V^coth(r'^)) (fT^tanh ij (C171) 

4/1 ^ 2 



A r'(<?) = - ^ (f7 - V?coth(r'V9)J (C172) 
This leads to the scaling form given in the text in equations (114) and (117). 

APPENDIX D: FINITE SIZE PROPERTIES 

In this Appendix we sketch the derivation of the finite size results from the finite size measure for the RFIM. 

If spins at both ends are fixed to have the same (resp. opposite) value, there are an odd (resp. even) number of 
bonds in the finite size measure. 

Assuming for definitness that spins at both extremities are fixed to the values So = +1 and Sl = — 1. Then domains 
closest to the boundaries are domain walls of type A which as Sinai walkers see reflecting boundary conditions. In 
this case, as explained in [28], the probability iVp fc ^ 2 (/i, Z 2 , ■ ■ ■ , hk+2) that the system at scale T has (2k + 2) bonds 
(k = 0, 1, . . .), with respective lengths (h, . . . , hk+2) is 

(2k+2 \ 
L - J2 M ( D1 ) 

where P^(l) are the bulk length distributions (68), 7r = Tp +l T , and the Ep(l) the length distribution of boundary 
bonds (see [28]). The normalisation is 
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]T / N^ 2 (l u l 2 ,...,l 2k+2 ) = l (D2) 

k = J l l'--- l 2k + 2 



The probability that the system has (2k + 2) domains at scale T is 

/2fe+2 
E+(l 1 )Pf(l 2 )...P+(l 2k+1 )Er(l 2k+2 )h5(L - ]T h) 
lf-hk+2 i = \ 

so that the Laplace transform with respect to the length of the generating function is 

E+{q/2g)E-{q/2g) 2g 1 - P+(q/2g)Pf(q/2g) 



>+(q/2g)Pf(q/2g) q 1 - zP+(q/2g)Pf (q/2g) 



(D3) 



(D4) 



where p = ^L. This leads to the results given in the text. 

The magnetization can also be obtained. In the "full" renormalized landscape, the magnetisation is given by 
m l = Ei=f +2 (-l) m ^- The probability that it has value M simply is 

fe=+oo 

F L {M-T)= / E + (h)p-(l 2 )P + (l 3 )..P+(l 2k - 1 )E-(l 2k+2 ) (D5) 

2fc+2 i=2k+2 

l r S(L - h)S(M - J2 ( D6 ) 

i=l i=l 

This leads to the result given in the text. 
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